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COUNTING COLOURED PLANAR MAPS: DIFFERENTIAL EQUATIONS 


OLIVIER BERNARDI AND MIREILLE BOUSQUET-MELOU 


Abstract. We address the enumeration of g-coloured planar maps counted by the number of 
edges and the number of monochromatic edges. We prove that the associated generating func¬ 
tion is differentially algebraic, that is, satisfies a non-trivial polynomial differential equation 
with respect to the edge variable. We give explicitly a differential system that characterizes 
this series. We then prove a similar result for planar triangulations, thus generalizing a result 
of Tutte dealing with their proper g-colourings. In statistical physics terms, we solve the 
g-state Potts model on random planar lattices. 

This work follows a first paper by the same authors, where the generating function was 
proved to be algebraic for certain values of q, including q = 1,2 and 3. It is known to be 
transcendental in general. In contrast, our differential system holds for an indeterminate q. 

For certain special cases of combinatorial interest (four colours; proper g-colourings; maps 
equipped with a spanning forest), we derive from this system, in the case of triangulations, an 
explicit differential equation of order 2 defining the generating function. For general planar 
maps, we also obtain a differential equation of order 3 for the four-colour case and for the 
self-dual Potts model. 


1. Introduction 

A planar map is a connected planar graph, given with one of its proper embeddings in the 
sphere, taken up to continuous deformation (Figure 1). The enumeration of planar maps is a 
combinatorial problem that has attracted a lot of interest since the sixties, in connection with 
graph theory [41, 26], algebra [28, 31[, theoretical physics [3, 16, 21[, and computational geom¬ 
etry [18, 37[. Several important approaches have been developed: recursive [40[, bijective [36[, 
using matrix integrals [16[, or using connections with the characters of the symmetric group [30[. 
We only give very few references, as a complete bibliography would take dozens of pages. 

From the combinatorial and physical point of view, it is natural to count planar maps equipped 
with an additional structure: for instance a spanning tree [34[, a proper colouring [45, 50[, an 
independent set of vertices [11, 13, 15[, a configuration of the Ising or Potts model [32, 23[, a 
self-avoiding walk [22[. (Again, we give very few of the relevant references.) The first result of 
this nature probably dates back to 1967 with Mullin’s enumeration of planar maps equipped with 
a spanning tree [34[. The second attempt is due to Tutte, who, in the early seventies, started 
to study maps — more precisely, triangulations — equipped with a proper colouring [45[. In 
the decade that followed, he devoted at least eight other papers to this problem [43, 42, 44, 46, 
47, 48, 49, 50[. His work culminated in 1982, when he proved that the series H{w) counting 
g-coloured rooted triangulations by vertices satisfies a polynomial differential equation [49, 50[ : 

2g2(l _ q)yj + [qw + 1077 - 6wH')H'' + - q)(20H - l8wH' + Qw^H”) = 0. (1) 
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We say that H{w) is differentially algebraic. Equivalently, the number hn of rooted triangulations 
with n vertices satisfies the following simple recurrence relation: 

n 

q{n+ l)(n + 2)/i„+2 = q{q - 4)(3n - l)(3n - 2)/i„+i + 2^i(i + l)(3n - 3i + l)hi+ihn+ 2 -i, 

i=l 

with the initial condition /i 2 = q{q — ^)- For instance, /13 = q{q—l){q — 2) is the number of proper 
g-colourings of a triangle. To date, this recursion remains entirely mysterious, and Tutte’s tour 
de force has remained isolated. 

Let us be more precise about the content of this tour de force. Tutte started from a generating 
function G{w; x, y) = G{x, y) which counts a larger family of g-coloured maps according to three 
parameters (as before, w counts vertices). The above series H{w) is G(r(;;l,0). Using the 
deletion/contraction properties of the chromatic polynomial, he easily established the following 
functional equation: 


2 , xy ^ 2 .^.G{x,y)-G{l,y) ^ ^G{x,y)-G{x,Q) 


G{x,y) = xq{q-l)w H- G{l,y)G{x,y)-x yw 

qw X — I 

Observe that the divided differences (or discrete derivatives) 


( 2 ) 


G{x,y)-G{l,y) G{x,y) - G{x,0) 

x-1 y 

prevent us from simply specializing a: to 1 and y to 0 to obtain an equation for H = G(1,0). 
In fact, the variables x and y, called nowadays catalytic variables [54, 11], are essential to write 
this equation, even though they are not very important combinatorially. The whole challenge is 
to determine the nature of the series H = G(1,0): is it algebraic, as the generating function of 
many classes of planar maps? Is it D-finite, that is, a solution of a linear differential equation, 
as the generating function of maps equipped with a spanning tree [34]? Is it at least differen¬ 
tially algebraic? Tutte answered the latter question positively by deriving from the functional 
equation (2) the differential equation (1) satisfied by H. 

It is known that the solutions of polynomial equations with one catalytic variable are system¬ 
atically algebraic [11[. Such equations commonly occur in the enumeration of families of planar 
maps with no additional structure. Moreover, in the past decade, progress has been made on 
linear equations with two catalytic variables, which one typically encounters when counting plane 
lattice walks confined to a quadrant. In this context, it is precisely understood when the generat¬ 
ing function is D-finite (see e.g. [7, 12, 33] and references therein). Hence the key difficulty with 
Tutte’s equation (2) is the occurrence of two catalytic variables, combined with non-linearity. 


A few years ago, we started to resurrect Tutte’s technique in order to solve a more general 
problem: the Potts model on planar maps. In combinatorial terms, this means that we count all 
g-colourings, not necessarily proper, but with a weight iz™, where v is an indeterminate and m 
is the number of monochromatic edges (edges whose endpoints have the same colour). The case 
iz = 0 is thus the problem solved by Tutte. Moreover, we do not only study degree-constrained 
maps (triangulations) as Tutte did, but also general planar maps (counted by edges and vertices). 

In a first paper on this topic [2], we wrote the counterpart of (2) for each of these two problems. 
Then, we performed a first step, by establishing an equation with only one catalytic variable, y. 
But this equation holds only for values of q of the form 2 -|- 2 cos(j 7 r/m), for integers j and m 
(with q 0,4). Moreover, the size of this equation grows with m. Nevertheless, equations with 
one catalytic variable are much better understood that those with two [ 11 ], and we were able 
to prove that the generating function of g-coloured maps is algebraic for all such values of g, 
including q = 2 and g = 3. It is known to be transcendental in general [2[. 

In this paper, we perform the second step, and prove that for g an indeterminate, the gener¬ 
ating function of g-coloured planar maps (or triangulations) is differentially algebraic. We give 
an explicit differential system that characterizes it. This system has the same form for general 
maps and triangulations. In some special cases (proper colourings; four colours; spanning forests; 
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self-dual Potts model) we derive from it an explicit differential equation of small order defining 
the generating function. In particular, we recover Tutte’s result (1) for properly coloured trian¬ 
gulations. For the convenience of the reader, all calculations are performed in an accompanying 
Maple session, available on the web pages of the authors. 

Here is an outline of the paper. In Section 2 we give definitions about maps, the Potts model 
and its combinatorial counterpart, the Tutte polynomial, as well as notation on power series. 
In Section 3 we state our main result and give an idea of Tutte’s sophisticated approach (some 
would say obscure) on a simple example. In Sections 4 and 5 we establish differential systems for 
g-coloured planar maps and g-coloured triangulations, respectively. We simplify these systems in 
Section 6. The next sections deal with special cases: proper colourings (Section 7), four colours 
(Section 8), then the q = 0 case, which corresponds to the enumeration of maps equipped with 
a spanning forest (Section 9), and finally the self-dual Potts model (Section 10). 

Our differential system for planar maps appeared (without proof) in a survey on map enu¬ 
meration published in 2011 [9]. 


2. Definitions and notation 


2.1. Planar maps 

A planar map is a proper embedding of a connected planar graph in the oriented sphere, 
considered up to orientation preserving homeomorphism. Loops and multiple edges are allowed 
(Figure 1). The faces of a map are the connected components of its complement. The numbers 
of vertices, edges and faces of a planar map M, denoted by v(M), e(M) and f(M), are related by 
Euler’s relation v(M) +i{M) = e{M) -|- 2. The degree of a vertex or face is the number of edges 
incident to it, counted with multiplicity. A corner is a sector delimited by two consecutive edges 
around a vertex; hence a vertex or face of degree k defines k corners. Alternatively, a corner can 
be described as an incidence between a vertex and a face. The dual of a map M, denoted M*, 
is the map obtained by placing a vertex of M* in each face of M and an edge of M* across each 
edge of M] see Figure 1. A triangulation is a map in which every face has degree 3. Duality 
transforms triangulations into cubic maps, that is, maps in which every vertex has degree 3. 

For counting purposes it is convenient to consider rooted maps. A map is rooted by choosing a 
corner, called the root-corner. The vertex and face that are incident at this corner are respectively 
the root-vertex and the root-face. The root-edge is the edge that follows the root-corner in 
counterclockwise order around the root-vertex. In figures, we indicate the rooting by an arrow 
pointing to the root-corner, and take the root-face as the infinite face (Figure 1). This explains 
why we often call the root-face the outer face and its degree the outer degree (denoted drf(M)). 

From now on, every map is planar and rooted. By convention, we include among rooted 
planar maps the atomic map having one vertex and no edge. 




Figure 1 . A rooted planar map and its dual (rooted at the dual corner). 
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2.2. The Potts model 

Let G be a graph with vertex set V{G) and edge set E{G). Let v be an indeterminate, and 
take g S N. A colouring of the vertices of G in g colours is a map c : V (G) — )■ {1,. .., g}. An edge 
of G is monochromatic if its endpoints share the same colour. Every loop is thus monochromatic. 
The number of monochromatic edges is denoted by m(c). The partition function of the Potts 
model on G counts colourings by the number of monochromatic edges: 

Pg(<Z,^)= E 

The Potts model is a classical magnetism model in statistical physics, which includes (when 
q = 2) the famous Ising model (with no magnetic field) [53, 52]. Of course, PG(q, 0) is the 
chromatic polynomial of G, which counts proper colourings (no monochromatic edge). 

It is not hard to see that Pg( 9) is a polynomial in q and v. We call it the Potts polynomial 
of G. Observe that it is a multiple of q. We will often consider q as an indeterminate, or evaluate 
Pciq, at real values q. 

We define the Potts generating function of planar maps by: 

M{y) = M{q, i., w,t-y) = -y] Pm(<?, (3) 

9 m 

where the sum runs over all planar maps M. Since there are finitely many maps with a given 
number of edges, and Pm((?, J^) is a multiple of q, the generating function M{y) is a power 
series in t with coefficients in Q[( 7 , v, w, y], the ring of polynomials in q, v, w and y with rational 
coefficients. The expansion of M at order 2 reads 

M{y) = w + {w‘^y‘^{q — 1 + v) + wyv) t + {2w^y^{q — 1 + + w^y^{q — 1 + i^^) 

+w^iy{y + 3y^){q - I + v) + wv‘^{y + y^)) t^ + O(t^), 

as illustrated in Figure 2. In combinatorial terms, M counts g-coloured planar maps by ver¬ 
tices, edges, monochromatic edges and outer degree, with the convention that the root-vertex is 
coloured in a prescribed colour (this accounts for the division by q). 



e = 0 

G' . O 

e = 2 


- O 

( 2 DC 


(1) 

(1) (iK*^ 


(4) v y 

PM{q) 

Q 

1 

{q-l + v) V 


{q-\W) 

v{q-l + u) 



Figure 2. The planar maps with e = 0,1, and 2 edges, and their Potts poly¬ 
nomials (divided by q). The maps are shown unrooted, and the numbers in 
parentheses indicate the number of corresponding rooted maps. 


2.3. Specializations 

We consider in this paper several specializations of the Potts polynomial. For some of them, 
the combinatorial meaning is obvious: for instance, Pg( 9 i 0 ) counts proper q-colourings, and 
Pg(4, p) counts 4-colourings by monochromatic edges. To understand the significance of some 
other specializations, it is useful to relate Pg to another invariant of graphs: the Tutte polynomial 
TG(/r, i^). It is defined as follows (see e.g. [4]): 

Tg(m,i^):= E (m-(4) 
SCE(G) 
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where the sum is over all spanning subgraphs of G (equivalently, over all subsets of edges) and 
v(.), e(.) and c(.) denote respectively the number of vertices, edges and connected components. 
The equivalence with the Potts polynomial was established by Fortuin and Kasteleyn [25] : 

PG(g,u)= ^ = TgM, (5) 

SCE(G) 

for q = {fj, — l)(u — 1). The Tutte polynomial satisfies an interesting duality property: if G 
and G* are dual connected planar graphs (that is, if G and G* can be embedded as dual planar 
maps) then 

=TG{v,fi). (6) 

This gives a duality relation on the Potts generating function defined by (3): 

M(( 7 ,u,w,t;l) = w^qM{q, n,{wq)~\tw{iy - (7) 

The connection (5) between Pq: and Tq allows us to understand combinatorially the limit q = 0 
of Pq /q: for a planar map M, 

lim - Pm(<7, v) = {v- T^f (1, v) 

9-j-o q 

^ (^_l)e(G) by (4), 

C connected on M 

= {v- l)fi^*)-i Tm* (u, 1) by (6), 

= (u-1)^^“*^"^ ^ by (4) again, 

F forest on M* 

where the first sum runs over all connected subgraphs C of M, and the second over all spanning 
forests F of M* (subsets of edges of M* with no cycle). Hence, if S{q,i',w,t) is the Potts 
generating function of some family of planar maps M (meaning that each map M is given a 
weight ^PM{q, and A4* denotes the set of duals of maps of M., 

5(0,l+/3,u;,t) = ^HG)^v(M)^e(M) 

MeM,C 

= \ ( 8 ) 

P M*eM-,F 

where the first sum runs over all maps M of Af equipped with a connected subgraph G, and the 
second over all maps M of J\4* equipped with a spanning forest F. 

Our final specialization is the self-dual Potts model, in which a map and its dual get the same 
weight. In view of the duality relation (7), this means that /r = u, so that q—{v— 1)^, and that 
w = 1/y^ = l/(u — 1). Denoting /3 = u — 1, we will thus consider 

+ 1, r\t- 1) = i ^ Tm(1 + /3,1 + /3) (9) 

^ M 

where the sum runs over all planar maps. One motivation for studying this specialization is 
that it should coincide with a critical random-cluster model [38, 19[, in which one chooses a 
map M and a subset S of its edges with probability proportional to q^l^) p- The 

partition function of this model is 

Y =9M(9,/3 + l,r'A;l). 

M SCE{M) 

The duality relation (7) implies 


M{q,/3,t)=q/3-^M{q,q/l3,t), 
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and one can thus expect the critical point (3c{q) to satisfy /3c(9)^ = 9- The critical partition 
function would then be given by (9). Note that the connection between criticality and self¬ 
duality has recently been proved for the regular square lattice [1]. 

2.4. Power series 

Let A be a commutative ring and x an indeterminate. We denote by A[x] (resp. ^[[a;]]) the ring 
of polynomials (resp. formal power series) in x with coefficients in A. The coefficient of a;" in a 
series F{x) is denoted by [x'^]F{x). If A is a field, then A{x) denotes the field of rational functions 
in X. This notation is generalized to polynomials, fractions and series in several indeterminates. 
For instance, the Potts generating function of planar maps M(q, z/, ic, t; y) defined by (3) belongs 
to Q[( 7 , w, y] [[t]]. To lighten notation, we often omit the dependence of our series in q, v and 
w, writing for instance M(t;y), or even M{y). 

If IK is a field, a power series F{t) G IK[[f]] is algebraic (over K(t)) if it satisfies a non-trivial 
polynomial equation P(t, F(t)) = 0 with coefficients in K. It is differentially algebraic if it satisfies 
a non-trivial polynomial differential equation P(t, P(t), F'{t), ..., = 0 with coefficients 

in K. 

For a series F in several variables x,y, ..., we denote by F^ the derivative of F with respect 
to x. 


3. Main result, and outline of the proof 

Our study of coloured maps started in [2] with the following equation with two catalytic 
variables x and y. 

M{x,y) = l+xywt{{v- l){y - 1) -f qy) M{x,y)M{l,y) 

+ xyffxv — l)M{x, y)M{x, 1) 

+ - ua.y) ^ 

X — 1 y ~1 

It defines uniquely a power series in t, denoted M(x,y) = M{q,v,w, t; x, y), which has polynomial 
coefficients in q,i/,w,x and y. The Potts generating function M{y) defined by (3) is essentially 
the specialization a: = 1 of M{x,y). More precisely, 

M{y) = M{q,v,w,t]y) = wM{q,iy,w,t; l,y) = wM{l,y). 

We then proved that when g ^ 0,4 is of the form 2 -|- 2cos{jTT/m), for integers j and m, the 
series M{y) satisfies an equation of the form: 

P{iy,w,t,y,M{y),Mi,... ,Mr) = 0, (10) 

for some polynomial P depending on j and m, and r (unknown) series in t, denoted by Mi = 
M(l),... ,Mr, which are independent of y [2, Cor. 10]. We call (10) a polynomial equation with 
one catalytic variable, y. In particular, when q = 1, every edge is monochromatic and (10) 
coincides with the standard functional equation obtained by deleting recursively the root-edge 
in planar maps [40]: 

M{y) = w + y'^tvM{yY + vty ^^^^'^ . (11) 

y- 1 

The classical way of solving (11) is Brown’s quadratic method, presented below in Section 3.1.1. It 
derives from (11) a polynomial equation satisfied by Mi. The quadratic method was generalized 
in ]11] to arbitrary (well-posed) polynomial equations with one catalytic variable: under minimal 
assumptions, their solutions are always algebraic. Using this approach, we proved in ]2] that 
for g = 2 -|- 2cos(j7r/m), the series M{q,v,w,t;y) is algebraic over M.{i',w,t,y) (and even over 
Q(g, p, w, t, y)). However, we could not obtain an explicit polynomial equation satisfied by M(jj), 
nor even by M(l) = Mi. We only constructed one for the three integer values of g of the above 
form, namely g = 1,2, and 3. And for g = 3, we were only able to do it when v = 0, that is. 
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when one counts proper 3-colourings. We expect the degree of Mi to grow with m. For q an 
indeterminate, the series Mi is transcendental. 

The main result of this paper is that when q is an indeterminate, and consequently for any 
real q, the series Mi is differentially algebraic over Q{q, i/,w,t). Moreover, we construct an 
explicit differential system that characterizes Mi. 


Theorem 1. Let q be an indeterminate, [3 = v — 1 and 

D{t, x) = (gu -I- - q{v + l)x + /3t{q - ‘i){wq + 13) + q. 


There exists a unique triple {P{t,x),Q{t,x),R(t,x)) of polynomials in x with coefficients in 
Q[( 3 ', u, w] [[f]], having degree 4,2 and 2 respectively in x, such that 


[x'^]P{t,x) = 1, P{0,x) = x'^(x — 1)'^, 

[x'^]R{t, x) = v + 1 — w{q + 2/3), Q(0, x) = x{x — 1), 


( 12 ) 


and 

Qm \pim) ~Rdx \ 1 mm) ■ 


(13) 


Let Pj{t) = Pj (resp. Qj, Rj) denote the coefficient of x^ in P(t,x) (resp. Q(t,x), R{t,x)). 
The Potts generating function of planar maps, Mi = M{q,v,w,t;l), can be expressed in terms 
of the series Pj and Qj using Mi := f^Mi and 


12 (/3^ -f qv) Ml -f P^ /4-|-2t(l-|-u — w(2/3 -I- q)) P 3 — P 2 2Qq = 4t (1 -1- w(3>j3 -\~ q)) ■ (14) 

An alternative characterization of Mi is in terms of the derivative of Mi: 

2 (/3^ -|- qv^ M'l -|- (1 -|- u — w (2/3 -I- q)) Pz/^ — Ri = 2 -|- 2Pw. 

The series Mi is differentially algebraic, that is, satisfies a non-trivial differential equation with 
respect to the edge variable t. The same holds for each series Pj, Qj and Rj. 

Eq. (13) looks like a partial differential equation in t and x, but it is in fact a system of nine 
differential equations in t written in a compact form. Indeed, its numerator reads: 

2Q'tPD - QPfD - 2QPD[ = 2i?;PD - RP^D - 2i?P£);. 

This is a polynomial in x of degree 8, and each of its nine coefficients gives a differential equation 
(with respect to the variable t) relating the eleven series Pj, Qj and Rj. For instance, extracting 
the coefficient of a;® gives : 

2Q^P4 - Q2Pi = 0. (15) 

The fact that we only obtain nine equations for eleven unknown series look alarming, but the 
initial conditions (12) give explicitly the two series P 4 and i? 2 - We also observe that (15), 
combined with the initial conditions (12), implies that Q 2 = 1. In fact, we will prove that the 
differential system (13), together with its initial conditions, determines the series Pj, Qj and Rj 
uniquely. For instance. 


Po = —-|- (q^w^ + IGfiw — Aqw -I- 8/3)t^ 

-I- 2{—j3q^w^ + q^w^ + 2j3qw^ — Aq^w'^ + l6/3^w + Af3qw + 2f3^ — 6 qw + Aj3)t^ A- 0{t‘^). 

Using (14), one can in principle construct a differential equation (DE) for Mi, but it would 
probably be very large. However, we will work out some special cases in details, and obtain for 
instance a DE of order 3 for four-coloured planar maps (Section 8). We predict the order of Mi 
to be 5 in general (Section 6.2). 

Our solution differs significantly from other published results on the Potts model, in that it 
ignores the catalytic variable y (which is set to its natural value 1) and describes the dependence 
of the series M(l) in the size variable t. In contrast, the results of [6, 29] give a precise description 
of the dependence of M{y) in y (in terms of elliptic functions), but the dependence in the 
size variable seems to remain elusive. This is this dependence that we have characterized, in 
differential terms. 
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3.1. A TOY example: uncoloured maps 

In this section, we illustrate on an example how one can derive differential equations, rather 
than polynomial equations, from equations with one catalytic variable like (10). Our example 
is Tutte’s equation for the enumeration of planar maps [40], counted by the number of edges 
(variable t) and the outer degree (variable y): 

M(t- y) = M{y) = 1 + ty^M{y f + ty ~ , (16) 

J/- 1 

where Mi stands for M(l) = M(t;l). This equation characterizes M{y) in the ring of formal 
power series in t. Indeed, the coefficient of t" can be computed by induction on n, yielding: 

M(t; y) = l + {y + y'^)t + (2y + 2y^ + Zy^ + 2y‘^)t^ + 0(f). (17) 

The reader can find this expansion, and much more, performed in the accompanying Maple 
session, available on the web pages of the authors. 

The standard method for solving (16) is Brown’s quadratic method (see [17], or [27, Sec¬ 
tion 2.9] for a modern account). It constructs from (16) a polynomial equation satisfied by Mi. 
We describe this solution in Section 3.1.1 below. We then present in Section 3.1.2 an alternative 
solution, which derives from (16) a system of differential equations satisfied by Mi. This gives 
the spirit of our construction of differential equations for g-coloured planar maps. The spirit, 
but not the details: in Section 3.2 we underline a few important differences between the solution 
of our toy example and that of g-coloured maps. 

3.1.1. The quadratic method: a polynomial equation for Mi. We first form in (16) a 
perfect square, as if we wanted to solve it for M{y): 

- i)^(y) + - y + 1 )^ = ( 2 / - i - y^^)^ - 4tj/^(y - 1 )^ + ^t^y'^iy - 1 )^ 1 - (I 8 ) 

Clearly, there exists a unique formal power series in t, denoted below Y, than cancels the left- 
hand side. Its nth coefficient can be computed by induction, using the expansion (17) of M(y): 

Y = l + t + 4ff + 2^1^ + 0{t^). 

Let us denote the right-hand side of (18) by 

^(y) = ^{t\y) = (y - 1 - y^if - ^ty'^iy - i)^ + 4t^y^(y- i)^i- 

This is a quartic polynomial in y, with coefficients in Q[[t]]. Since Y cancels the left-hand side 
of (18), it also cancels its right-hand side. Thus T is a root of A(j/). By differentiating (18) with 
respect to y, we see that Y is also a root of Ay, and hence a double root of A. 

Since A{y) has a double root, its discriminant (with respect to y) vanishes. This gives a 
polynomial equation satisfied by Mi: 

27ffMf + {1 - 18t) Ml + 16t - 1 = 0. 

Of course this can be easily solved, and we obtain the well-known generating function of planar 
maps counted by edges [40] : 

(1 - 12f)3/2 _ 1 + I 8 t 

3.1.2. An alternative approach: a differential system. We now describe another conse¬ 
quence of the fact that A(y) has a double root, this time in terms of differential equations. The 
procedure looks cumbersome, compared to the neat fact that the discriminant of A{y) is zero. 
But, applied to coloured planar maps, it will yield a system of differential equations that looks 
much nicer than any polynomial equations we could possibly compute. 

Denote Z = 1/T. Since T is a double root of A{y), there exists a polynomial P{t\y) = P{y) 
of degree 2 such that 


A(y) = (1 - yZfPiy). 
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Since Z is a formal power series in t, the same holds for the coefficients of P{y). Similarly, there 
exists a polynomial Q{y) of degree 2 such that 

^'yiy) = i^-yZ)Qiy), (19) 

and a polynomial R{y) of degree 3 such that 

A[{y) = {l-yZ)R{y). (20) 


We now want to eliminate A and Z in these three equations so as to find a (differential) relation 
between P, Q and R. We first eliminate Z, writing 


= — and —— 

A P A 


”P ' 


( 21 ) 


We now differentiate the first equation with respect to t, and the second one with respect to y: 


2A;A" ,A - A'/A( d fQ^ 
A2 dt\P 


2A(A" ,A - A(^A; _ d 
A2 dy 


~P 


The ratio of the left-hand sides is A'/A(, which, according to (19-20), is also Q/R. This gives: 


1 d [Q 


Q dt \ P 


R'^ 


~P 


Note the striking analogy with our differential system (13) for coloured maps. As in the coloured 
case, this equation gives, in a compact form, a system of differential equations in t relating the 
coefficients of P{y), Q{y) and R{y). We will not discuss which additional properties are needed 
to characterize them uniquely, but let us express the series Mi in terms of them. The first 
equation of (21) reads 

pa;" = q^a. 

This is a polynomial in y of degree 8. Extracting the coefficient of i/® gives the identity 

64^2P 2 M 1 -f let^Pa -Ql- 64tP2 = 0, 


which determines Mi in terms of the coefficients Pj and Qj of P{y) and Q{y). 


3.2. Some differences between the toy example and coloured maps 

Before embarking on the proof of Theorem 1, we want to underline a number of differences 
between the treatment just applied to uncoloured maps and our treatment of coloured maps 
below. Let us list the main steps of our approach. 

• We take q = 2 2cos(27r/m), and start from a polynomial equation in one catalytic 

variable satisfied by M{y), obtained in [2]. This is (25). 

• This equation has some similarities with (18): in the latter equation, the right-hand side 
is a polynomial in y with coefficients in Q[[t]], while in (25), the right-hand side is a 
polynomial in M{y) (or rather, in a series /(y) which contains the same information as 
M{y)), with coefficients in Q(( 7 , u, w)[[t]]. Hence the roles of y and M{y) are exchanged. 

• We derive from (25) that a certain polynomial A(x), of degree 2 m and involving m — 2 
unknown series in f, has m — 2 double roots Ii,, Im- 2 - This is Lemma 4. 

• We derive our differential system from this property. The counterparts of (19) and (20) 
do not have left-hand sides A; and A(, but variants of these two polynomials (Proposi¬ 
tion 5). 
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4. Differential system for coloured planar maps 

4.1. An equation with one catalytic variable 

Our starting point is an equation of [2]. We take q = 2 + 2cos(2fc7r/m), with k and m coprime 
and 0 < 2k < m. We write /3 = p — 1. We introduce the following notation: 

• I{t,y) = I{q,v,w,t;y) is a variant of the generating function M(y) = M{q,iy,w,t;y) of 
g-coloured planar maps defined by (3): 

I{t,y)= tyqM{y) + - - - + . (22) 

y y-1 

• N{y,x) and D{t,x) are the following (Laurent) polynomials: 

N{y,x) = l3{A-q){y- 1) + {q + 2/3)x- q, (23) 

with y = 1 /y, and 

D{t, x) = {qiy + P‘^)x'^ - q{iy + l)x + /3t(g - 4)(wg + /3) + g. (24) 

Let Tjn be the mth Chebyshev polynomial of the first kind, defined by 

Tm(cos0) = cos{m9). 

Then there exists to + 1 formal power series in t with coefficients in Q{q, v, w), denoted C'o(t),..., 
Cm{t), such that 

D{t,I{t,y)r/^ (25) 

\2y/D{t,I{t,y)) J ^ 

This is a combination of Corollary 10 and Lemma 16 from [2]. We call (25) the invariant 
equation, since it is derived from a certain theorem of invariants, in Tutte’s terminology [51]. 
We find convenient to denote 

771 

C{t,x) =J2^r{t)x^. (26) 

r—0 

Since Tm(n) is even (resp. odd) in u if m is even (resp. odd), the left-hand side of (25) only 
involves integer powers of D. 

As the invariant equation may look intimidating, let us consider an example. 

Example: bipartite maps. Let us take q = 2 (that is, to = 4 and k = 1), w = 1 and 
p = 0. Then M{y) is simply the generating function of bipartite maps, counted by the edge 
number (variable t) and the outer degree (variable y) . Let us show that in this case, the invariant 
equation (25) coincides with the standard equation obtained by deleting the root-edge. 

With our choice of q, v and w, we have 

fV(?/, x) = —2y, D{t, x) = x^ — 2x 2t + 2, 

while Tjn{u) = T 4 (m) = 1 — 8u^ + 8u^. The invariant equation thus reads 

4 

{I{yf - 21{y) + 2t + 2f- 8f{I{yf - 21{y) + 2t + 2) + Sf = I{yY ■ 

r—0 

In this equation, let us replace I{y) by its expression (22) in terms of M{y) (with q = 2), 
and then expand the equation in the neighbourhood of y = 1. Saying that the coefficients of 
{y — 1)“"^,..., (y — 1)° must vanish gives the values of the series Cr- 

Ci = 1, Ca = -4, Ca = At, Ci = 8(1 -k t) (27) 

and 

Co = -4:-A0t-Af + 32YM{1). 

In the invariant equation, let us replace each series Cr by its expression: we obtain 
yhiy^ - l)M{yf + {1 - y^ + yh)M{y) - ty^M{l) + y'^-I = Q. 
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Equivalently, 


M{y) = 1 + + ty'^ 


M{y) - M(l) 

- 1 


which can be obtained by deleting the root-edge in a bipartite map [40]. 


□ 


Let us return to the general case q = 2-|-2cos(2fc7r/m). By expanding (25) near y = 1, we can 
express the series Cr{t) in terms of the derivatives of M with respect to y, evaluated at y = 1. 
We will need the expressions of Cm, ■ ■ ■, Cm- 3 - Note that the first three are completely explicit 
in terms of q, v, t, w and m. The fourth one involves the unknown series M{1). 

Lemma 2. Denote P = v — 1, 5 = fP + qv and 7 = . We have: 

Cm = 5^'^ T^( 7 ), 

Cm-i = ( 2 m(u + 1 )V^ T,„( 7 ) + fi{q - 4) T:„( 7 )) , 

Cm -2 = ^ (2maoV^ T^( 7 ) + g/3(<z - 4)ai T:„( 7 )) , 

Cm -3 = { 2 mboVS Tmil) + fi{q - 4)61 T:„(7)) + 

with 

ao = 2tpS{q-4){l3 + wq) + q{m-\) ((l + v^) q - 2/3^) , 

Oi = 2tS {ly + 1 — w {2 P + q)) + q {m — 1) (u -|- 1) , 

60 = 6tw (q — 4) P 6 ((/3^ — g) TO -I- g (u -I- 1)) — 6 P^tS (to — 1) (u -I- 1) (g — 4) 

- g (to - 1) (to - 2) (u -I- 1) ((u^ - u -I- 1 ) g - 3/3^) , 

bi = 6twS (^q{m — l)(u^ -|- 2 u -|- g — 3) — (2 /3 -|- g) (5) 

— 6t(5 (to — 1) (g(l -I- — 2/3^) -|- q{m — 1)(to — 2) (/3^ — + v + l)g) ■ 

Example: bipartite maps (continued). When k — 1, m = 4, q = 2, w = l and u = 0, we 
have P = —1, 5 = 1 and 7 = 0 . Moreover, T,„( 7 ) = 1 and ’T'mi'l) — We can then check that 
the above lemma specializes to (27). In particular, Ci = Cm -3 does not involve M(l) because 

TL(7) = 0. □ 

Proof of Lemma 2. We multiply the invariant equation (25) by {y — 1)"*, so that it becomes 
non-singular at y = 1, and expand it around y = 1. We extract successively the coefficients of 
{y — 1 )°,... ,{y — 1 )^ and obtain in this way expressions of Cm, ■ ■ ■, Cm- 3 - For instance, since 

(y- l)Iit,y) =t + 0{y- 1 ), 

(y - ^)N{y,I{t, y)) = (q + 2P)t + 0(y - 1), 

{y - lfD{t,I{t, y)) = Sf + 0{y - 1), 
extracting the coefficient of {y — 1 )° gives 

b^lH^Tm(n)=Cmr, 

from which the expression of Cm follows. As we extract the coefficients of higher powers of 
(y — 1), derivatives of T^, taken at the point 7 , occur. We systematically express them in terms 
of Tm ( 7 ) and T'm ( 7 ) using the differential equation 

(1 - u^)T" (n) -uT'^{u) + to^T„(u) = 0. 

Similarly, derivatives of (1 — y)I{t,y) with respect to y occur, and this is why the expression of 
Cm -3 involves the series M(l). a 
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4.2. Some special values of y 

The invariant equation (25) reads 

Ratit,y,I{y),Co,...,C^)=0 (28) 

where Rat is the following rational function, with coefficients in Q{q, v,w): 

Rat(t, 2 /,a;,Co,...,c^) = D(f, x)™/^ 

\WD(t,x)J ^^0 

Observe that the first term of Rat is the only one that involves y. In our toy model (Section 3.1), 
we were dealing with (18), which read Rat(f, y, M{y), Mi) = 0 with 

Rat{t,y,x,c) = {2ty^{y- l)x + ty'^ -y+l^ - {y-I- y^t)'^ + Aty'^{y - 1)^ - At‘^y^{y- l)c. 

There, the first term was the only one that depended on x. To solve this toy equation, we 
considered a series Y = Y(t) cancelling this term, or equivalently, satisfying 

Ra4(t,y,M(y),Mi) = 0, 

where the derivative is taken with respect to the third variable of Rat. Similarly, we are going 
to solve (28) by cancelling Rat^ (this change from x to y explains why we wrote in Section 3.2 
that the roles of y and M{y) are exchanged). More precisely, we focus on 

Denoting n = [to/ 2J, the set of roots of T(^ is easily seen to be: 

{cos(j7r/m),0 < j < m} = {±cos{jTr/m),j = 1,... ,n} . (29) 

Hence we are interested in the equation 

N(Y,Iit,Y))^ = 4cosUTT/m)^Dit,I{t,Y)), (30) 

for some j G |1, n]. 

Lemma 3. Let us assume k = 1, so that q = 2 + 2cos(2'!T/m). There exist m — 2 distinct formal 
power series in t, denoted Yi,... that have coefficients in K(p, tc), constant term 1, and 

satisfy (30) for some j G |l,n]. 

Let us denote Ii(t) := I(t,Yi). This is a formal power series in t with coefficients in 
The m — 2 series Ii{f) are distinct, and for 1 < i < m — 2, 

dl 

L(0)^{0,1}, D{t,L)^0, —{t,Y,)^0. 

Example: one-coloured maps. We take <7=1, that is, m = 3 (and n = 1). The only relevant 
value of J is 1, and we thus want to solve 

N{Y,I{t,Y))^=D{t,I{t,Y)), 

for Y a series in t with constant term 1. We write Y = 1 + tZ. Using the definitions (23) and (24) 
of N and D, the equation satisfied by the series Z can be written as 

Z = V + t Pol {v,w,t, Z,M{1 +tZ)), 

for some polynomial Pol with integer coefficients. This equation allows us to compute by induc¬ 
tion on n the coefficient of in Z, and shows the existence and uniqueness of Z (and Y). More 
precisely, 

Y = l-\-vt + 2v^{w + \)t^ + p^(5 -I- lAw + 6w^)t^ + 0{t‘^). 
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This gives 


Ii{t) := I{t,Y) = 1/u + - w -l)t + 

i?(t,/i) = (l-l/u)" + 0(i), 
dl, , 1 ^, 

so that the last three properties of the lemma hold. 


□ 


Proof of Lemma 3. We denote qj = 4 cosQ'tt/to)^ = 2 + 2 cos{2jTT/m). Observe that q = qi- We 
assume that m = 2n + 1 is odd. The proof is similar in the even case. 


Existence of the series Yi. We are looking for series Y = Y{t) solutions of (30), of the form 
Y = 1 + tZ for some series Z = Z{t). Multiplying (30) by Z"^ gives the equation ^j{t, Z) = 0, 
where 

$j(t, z) = z^N{l + 1 + tz))'^ — qjZ^D{t,I{t, 1 + tz)). 

We are interested in non-zero solutions (since for Z = 0, we have Y = 1 and I{t,Y) is not well 
defined). Observe that zl{t, 1 + zt) is a formal power series in t with coefficients in K[u,?c,z]. 
Hence the same holds for zN{l + tz, + tz)) and z'^Dft, I(t,l + tz)), and finally for $j(t,z). 
Recall that u = /3 + 1. Expanding I, N and D at first order in t gives 


zl{t,l + tz) = 1 + 0{t), (31) 

zN{l + tz,I{t, 1 + tz)) = q + 2j3 — qz + 0{t), 

z‘^D{t,I{t, 1 + tz)) = ql3 + q + - q{f3 + 2)z + qz^ + Oft). (32) 

The equation ^jft, Z) =0 thus reads 

Pj{Z)+tSj{t,Z)=0, (33) 

where 

Pj{z) = (g + 2/3 - qz)^ - qj (g/3 + g + /3^ - g(/3 + 2)z + qz^) 
and Sjft,z) is a power series in t with coefficients in w, z]. By (33), the coefficient of 
in Z, denoted zq, must satisfy Pj{zo) = 0. Equivalently, zq = 1 + Pv with 

q{q - qj)v‘^ + q{qj - 4)u + 4 - q, = 0. (34) 


This equation has degree 1 in u if = g, that is if j = 1, and degree 2 otherwise. Given the 
expressions of qj and g = gi, its roots are found to be 




sin(j7r/m) 


and 


sin()7r/m) 


V ■ =■ - 

^ 2 cosfK/m) sin((j — 1 ) 71 / 771 ) ’ 


2 cosfK/m) sin((j + 1 ) 71 / 777 ) 

for / ^ 1. For / = 1 the only root is vff = 1/g. The roots v)j~ and v~ are distinct (we use here 
the assumption that m is odd). 

Now having fixed zq = 1 + Pv with v = v^, let us return to (33) and extract from it the 
coefficient of for p > 1. This gives 


Pj{zo)[f^\Z + (expression involving only \f]Z for i < p) = 0. 


Since v is not a double root of Pj , this allows us to compute the coefficient of fP in Z by induction 
on p, and thus determines Z completely. Moreover, \tP]Z has a rational expression in v and w 
(with real coefficients). 

It is easy to see that 


0 < u/ < • • • < 


1 


V' = 


2cos(7r/777) 


< 7;„ < ■■ ■ < V. 


2 : 


so that we have exactly {m — 2) distinct values . Let us denote them vi,, Vm- 2 - They give 
rise to exactly fm — 2) distinct series Y satisfying E(0) = 1 and (30) for some j € |1,77]. We 
denote them W,..., Em- 2 ) with = 1 + (1 + Pvi)t + 0{t^). We note that Vi is a real number. 
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while j3 = V — \ where v is, & formal parameter. Hence our m — 2 series Yi all differ from the 
constant 1, and I{t, Yi) is well-defined. We have thus proved the first statement of the lemma. 


Properties of Yi and Jj I{t,Yi). Let us prove that the series li := I{t,Yi) are distinct. 
Returning to (31) gives li = {1 + I3vi)~^ + 0{t). Hence the (m — 2) constant terms 1^(0) are 
distinct since the (m—2) numbers vi,..., Vm -2 are distinct. Moreover 1^(0) is obviously non-zero, 
and it is distinct from 1 because Vi ^ 0. 

Let us now prove that D{t,Ii) ^ 0. Returning to (32), we find 

{1 +l3vi)‘^D{tJi) = 0^{qv'^i - qv^Yl)+ 0{t) = +0{t), (35) 

V sm((j ± l)7r/m)/ 


where the value of the denominator depends on whether Vi is of the form v'^ or v~. At any rate, 
this is non-zero. 

We finally check that 


dy 


{t,Yi) 


1 

t{l + /3Ui)2 


+ 0 ( 1 ) 


is non-zero, and this completes the proof of the lemma. 


4.3. Some polynomials with common roots 

We still assume q = 2 + 2 cos(27r/m). Recall that D{t,x) and C{t,x) are polynomials in x 
with coefficients in Q((jf, p, i(;)[[t]], given by (24) and (26) respectively. Their degrees in x are 2 
and m, respectively. The coefficients of D are explicit, but those of C are unknown, apart from 
the three leading ones (Lemma 2). We now prove that several polynomials, related to C{t,x) 
and D(t, x), admit the series A of Lemma 3 as common roots. 


Lemma 4. Each of the series li := I(t,Yi) defined by Lemma 3 satisfies 


C{t,h)^ = D{t,L)^, (36) 

m f)D 

D{t,L)^it,Ii) = yC(t,L)|^(LA), (37) 

m f)D 

D{t,L) — it,L) = -C{t,I,) — {t,h). (38) 


Since D{t,Ii) 0, it follows from (36) and (37) that L is actually a double root of — D™. 


Proof. Let us denote 


G{t,y,x) 


N{y,x) 

2^D{t,x) 


Recall that D{t, li) is a formal power series in t with coefficients in K(p, w). By (35), its constant 
term is of the form /3'^r^/(l+Pv)^, for some real numbers r ^ 0 and v. So we can define •\/ D{t, li) 
as a formal power series in t with coefficients in K(p, w). 

By construction of the series Yi, we have G{t, y, li) = ± cos{j'K/m), for some j G |1, nj. Given 
that Tm(cosa;) = cos{mx), each series Yi satisfies 


T„ {G{t, Y, I{t, Yi))) = e and {G{t, Y, I{t, Yi))) = 0, 


where e = ±1. Hence the invariant equation (25) gives 


G{t,Ii)=eD{t,Ii)'^G 


(39) 


from which (36) follows. 
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Let us now differentiate the invariant equation (25) with respect to y: 

/fin fin f)T \ 

+ D{t,I{t,y))'^/^ \^ — {t,y,I{t,y)) + —{t,y,I{t,y)) — {t,y)j {G{t,y, I{t,y))) 

dG, , . 

or, if we want to keep things in a reasonable volume, 

T,„(G) + (g; + GU;) TUG) = G'J'y 

When y = Yi, then Tm(G) = e, T'UG) = 0 and ly 0 (by Lemma 3), so that 

m 


= g;. 

Combined with (39), this gives (37). 

Similarly, differentiating the invariant equation with respect to t gives 

^ ^ j^ra/2 (g, ^ ^ ^ 

When y = Yi, then Tm(G) = e, T'UG) = 0, and G(, is given by (40). Hence 

= GL 


(40) 


Combined with (39), this gives (38). 


Proposition 5. Let L := I(t,Yi) be the (m — 2) series defined in Lemma 3. There exists a 
triple (P(t, a:), Q(t, cc), i?(t, cc)) of polynomials in x with coefficients in K(u, r(;)[[t]], having degree 
at most 4, 2 and 2 respectively in x, such that 

m-2 

Git, a;)^ - Dit, x)"" = Pft, x) (a; - lU, (41) 

m-2 

Dit,x)G'U^x)-^D'U,x)Cit,x)=Qit,x) ]J(x-/i), (42) 

m—2 

Dit,x)Gi{t,x) - ^Di{t,x)G{t,x) = R{t^x) J|(x-/i)- (43) 

i=l 

Proof. The polynomial G(t, x)^ — 0(1, x)™ has degree (at most) 2m in x, and admits each of the 
(to — 2) distinct series L as a double root (Lemma 4). Thus, there exists a polynomial P{t,x) 
of degree at most 4 such that (41) holds. Let us write 

m—2 ^ 

Pit, x) = {Git, x)2 - Dit, x)U p(i-xrU ' 

Recall that G and D are polynomials in x with coefficients in K(u, w)))^]], that the series L 
belong to K(u, r(;)[[t]] and have a non-zero constant term (Lemma 3). Expanding in x the above 
expression thus proves that the coefficients of P also belong to K(u, r(;)[[t]]. 

Similarly, Il(f, x)G(,(f, x) — ^I?(,(t, x)G(t, x) is a polynomial of degree at most to in x, since 
£)(<, x)G(.(<, x) and yZ?(,(f, x)G(t, x) are polynomials of degree to -|- 1 having the same leading 
coefficient. Hence by (37), there exists a polynomial Qit, x) in IR(u, ■u;)[[t]][x], of degree at most 2, 
such that (42) holds. 

Finally, given that D[it,x) has degree 0 in x, and G((t,x) degree at most m — 2 (Lemma 2), 
the polynomial Zl(t,x)G((t,x) — x)G(t,x) has degree at most to. Equation (38) implies 
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the existence of a polynomial R{t,x) in K(i/, of degree at most 2, such that (43) holds. 


4.4. Differential system 

We still assume that q = 2 + 2cos(27r/m). We will prove that the series P, Q and R of 
Proposition 5, suitably normalized into series P, Q and R, satisfy all equations of Theorem 1. 
We repeat this theorem for convenience. 


Theorem 1 (repeated). Let q be an indeterminate, P = v — 1 and 

D{t, x) = {qv + P‘^)x^ - q{v + l)x + pt{q - A)(wq + P) + q. (45) 


There exists a unique triple {P{t,x),Q{t,x),R(t,x)) of polynomials in x with coefficients in 
Q[( 7 , F, w] [[t]], having degree 4,2 and 2 respectively in x, such that 


and 


[x^]P{t,x) = 1, 

[x‘^]R{t, x) = + 1 — w{q + 2/3), 


1 d f \ _ 1 d 
Qdt \¥^) ~ Rdx 


P{0,x) 

= x'^{x — 1)2, 

(46) 

Q{0,x) 

= x{x — 1), 

/ Ft^ \ 

1 PP2 ) 


(47) 


Let Pj{f) = Pj (resp. Qj, Rj) denote the coefficient of x^ in P(t,x) (resp. Q(t,x), R{t,x)). 
The Potts generating function of planar maps, Mi = M(q,i',t,w,l), can be expressed in terms 
of the series Pj and Qj using Mi := f^Mi and 


12 {P^ + qv) Ml + P|/4 + 2t (1 + F - wifip + q)) P^-P 2 + 2Qo = 4t (1 + u>(3/3 + q)). (48) 

An alternative characterization of Mi is in terms of the derivative of Mi: 

2 {P^ + qv)M'i + {l + v-w (2/3 + q)) P 3/2 - i?i = 2 + 2Pw. (49) 

The series Mi is differentially algebraic, that is, satisfies a non-trivial differential equation with 
respect to the edge variable t. The same holds for each series Pj, Qj and Rj. 


A differential equation relating P, Q and R. We first prove that (47) holds for P, Q and 
R. This results from the elimination of C and ni(^ ~ ^i) ^^e three equations of Proposition 5. 

The derivation is analogous to that of Section 3.1.2. Let us first eliminate the product, in such 
a way that Q‘^/{PD'^) and Rf /[PD'^) naturally appear: 

i:)2(c'2 _£,m) - -p^ -ffyffcTZTjffLy ~ pzp- 

Let us differentiate the first equation with respect to t, and the second one with respect to x. 
The ratio of the two resulting identities is 

pc; - f 
DC[ - fD[C ’ 

which, according to the last two equations of Proposition 5, is Q/R. This gives 


]_d_(\ 

Qdt ypD^) ~ Rdx ypD'^) ' 


(50) 


This equation coincides with (47) but with hats over the letters P, Q and R. The series P, Q 
and R occurring in the theorem will simply be normalizations of P, Q and R by multiplicative 
constants independent from t and x. 
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The leading coefficients of P, Q, R. Proposition 5 gives 
[x^]P{t,x) = {C{t,xf - , 

[x^]Q{t,x) = [x'^] [D{t,x)C'^{t,x) - ^D'^{t,x)C{t,x)^ , 
[x‘^]R(t,x) = [x™] (^D{t,x)C't(t,x) - ^D[{t,x)C{t,x)^ . 


Recall that D is defined by (45), and that Lemma 2 gives the leading coefficients of C{t, x). This 
allows to determine the leading coefficients of P, Q, R: 

rP 4 [x^]P{t,x) = {Tmilf - 1) = ;^( 9 / 4 - 1)(5— 1 t:„( 7 ) 2 , 

I 02 := [x^]Q{t, x) = q(3{q/A - l) 5 ^rn-i )/2 ( 51 ) 

[r 2 := [x'^]R{t,x) = qPiqjA- T(„( 7 )(u + 1 - w{q + 2P)), 


where 6 = 13^ + qv and 7 = (g + 2/3)/ {2VS). In the expression of P 4 , we have used the fact that, 
for ctll X 

{u^-l)TUuf=m^TM^-l). (52) 

Observe that these coefficients are independent of t. Let us define 


Pit,x) = 


Pit,x) 


Qit,x) = 


Qit,x) 


and R{t,x) = 


R{t, x) 


(53) 


P 4 Q 2 Q 2 

We have intentionally normalized R{t,x) by Q 2 rather than R 2 . It then follows from (50) that 
P, Q, R satisfy the differential system (47). Moreover 

[a;^]P(t, a;) = 1, [x^]Q{t,x) = l and [x'^]R{t,x) = v + 1 — w{q2^), (54) 

so that the left-hand side of the initial conditions (46) hold. (We do not give the value of 
[x^]Q{t, x) in the statement of the theorem because it is a consequence of the differential system 
and the initial conditions, as will be seen in Section 4.6.) 


The case t — Q. Let us now establish the right-hand side of the initial conditions (46) by 
determining the values P{0,x) and Q(0,a;). Proposition 5 gives 


C{0,xf - D{0,x)^ 

T71 

D{0,x)C',i0,x) - -D'^{0,x)Ci0,x) 


m—2 

P{0,x)l[{x-I,{0)f, 

m—2 

Q{0,x) {x- li{0)). 

2=1 


(55) 

(56) 


We will now combine our knowledge of C{Q, x) obtained in [2] with the properties of the series 
Ii{t) gathered in Lemma 3 to determine P{Q,x) and Q(0,x). 

Recall that Cr{t) denotes the coefficient of x'^ in C{t,x) (see (26)). The constant term of 
the series Cr{t) has been determined in [2, Lemma 16]. With the notation used there, C'r(O) = 
Lr( 0 ; 1 ), where 


f2 Lr{t-, y)x^ = P(f, xT/^ 

7^0 \^\/D{t,x) 


where N and D are given by (23) and (24). Hence 


where N = N{1, x) 
first want to factor 


m / \r \ 

C{0,x) ^ J2CriO)x'^ = 

(g -|- 2/3)x — q and T) = P(0, x) = {qv + I5'^)x^ 


(57) 

q[v -I- l)a; + q. Thus we 



C{{),xf - D{Q,x)'^ = V'^ 
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Denote n = [to/ 2J. By (52) and (29), 


1 / ■ \2 
Tm{uf -1= - n (m"-cos2^) 

i=i ^ ^ 


^ ■ ( 58 ) 


(We have also used the fact that the dominant coefficient of Tm is 2™ ^.) Thus 

1 


^(O,^)^ -D(0,x)™ = -{Af^ - AV){Af^)^'" 


= ^(g-4)(x-l)2(AfY^ 


L^J . . X 2 

n (Af2-4I?cos2^ 

/ . X 2 

]J (aA'^-AVcos^ — 

j=i ^ ^ 


(59) 


Let us now compare this to (55). By Lemma 3, li{0) ^ 1 for all i, so that (x — 1)^ is necessarily 
a factor of P(0, x). We will now prove that x^ is also a factor of P(0, x). 

Consider the term obtained for j = 1 in (59). Using 4cos^(7r/TO) = q, we get 

— 4Pcos^ — = — qV = xj3{q — 4)(g — xiq + j3)), 

m 

which has a factor x. By Lemma 3, /i(0) ^ 0 for all i, so that x^ is necessarily a factor of P(0, x). 
We have proved that P(0,x), hence also P{0,x), is divisible by x^{x — 1)^. Moreover P{0,x) 
has degree 4 and constant term 1, thus P(0,a:) = x'^{x — 1)^. 

We now wish to determine Q{0,x). We have just seen that a: = 0 and x = 1 are dou¬ 
ble roots of (7(0, a;)^ — P(0,a;)'". So they cancel as well the derivative 2(7(0, a;)(7(0, a;)^ — 
mD{0, a:)'"“^P(,(0, x), and since P(0,0) 0 and P(0,1) ^ 0, they must also cancel P(0, x)C'^{0,x) 

^P(,(0,a;)(7(0,a;). Let us now return to the factorisation (56). Since Ii{Q) yf 0,1, we see that 
x{x — 1) must divide Q{0, x). Hence it also divides (3(0, x). Moreover we have proved above that 
[x’^]Q{0,x) = 1 so that Q{0,x) = x{x — 1) as stated in the theorem. 


4.5. The Potts generating function of planar maps 

We still assume that q = 2 + 2cos(27r/m). Let us prove that the Potts generating function 
Ml is related to the Pj’s and Qj’s by (48) and (49). It follows from the first two identities of 
Proposition 5 that 

P (pc; - f (c2 _ D^) , 

where all series and polynomials are evaluated at {t,x). Using the normalization (53), this gives 

p (PC; - |p;c) ' = qm^/3^(q/4 - 1)Q^ {C^ - D^) , (60) 

since Q\ = qm^I3'^{q/A — l)Pi by (51). Recall that D is defined by (45), that the leading 
coefficients of (7 are given in Lemma 2, and that (7^ — P™ has degree 2m. Recall also the known 
values (54) of P4 and Q 2 . Extracting from (60) the coefficients of and 

gives: 

• for the coefficient of a tautology, equivalent to (52) taken at u = 7, 

• for the coefficient of an interesting relation between P 3 and Qi, namely 

P3 = 2Qi 4f(l -f p) — 4:tw(2/3 q), (61) 

• for the coefficient of the identity (48). In the calculation we use once again (52) 

to relate Tmij) and as well as (61) to express Qi in terms of P3. 
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The second characterization (49) of Mi is obtained in a similar fashion by combining instead 
the second and third identities of Proposition 5: they imply 

R = Q [dC', - jD',c) , 

and extracting the coefficient of gives an expression of M( in terms of Qi and i?i, which 

we transform into (49) using (61). 

4.6. Uniqueness of the solution 

The arguments in this subsection apply whether q is an indeterminate, or q = 2 + 2 cos(27r/m). 
The differential system of Theorem 1 can be written as 

2Q'tPD - QP^D - 2QPD[ = 2R'^PD - RP'^D - 2RPD'^. (62) 

Since P, Q and R have respective degree 4, 2 and 2 in x, this identity relates two polynomials in 
X of degree at most 8. Recall that P 4 = i ?2 = 1- Extracting the coefficient of x® gives Q2{t) = 0, 
which, with the initial condition Q2(0) = 1, implies Q 2 {t) = 1. Hence the leading coefficients 
of P, Q and R (that is, the series P 4 , Q 2 and R 2 ) are independent of t, and the left-hand side 
of (62), as well as its right-hand side, has degree at most 7 in x. And we are left with eight 
unknown series. 

We denote Pij := [f^]Pj, and similarly for Q and R, so that 

P{t, x) = PijfxP 

Let Ci be the following 8-tuple of coefficients: 

5 Qi,0^ Qi,l: Ri—1,0: Rz—1,1) ■ 

The right-hand side of (46) gives us the values of P{t,x) and Q{t,x) at t = 0, so that 

Co = (0,0,1,-2; 0,-1; 0,0). 

We will show by induction on z > 1 that the differential system determines the eight coefficients 
of Ci, and that these coefficients are rational functions of q,/3,w. 

For z > 1 and 0 < j < 7, the equation Eq^ ^ obtained by extracting the coefficient of C~^x^ 
in (62) reads 

(2*1 ~ *2 ~ 2z3)(5ij J3 “ (2jl ~ J2 ~ 

•*1+^2+ ^3='* il+i2+i3=i-l 

Jl+72+73=.7 Jl+j2+J3=J + l 

where Di j = \Cx^]D[t,x). This is a linear equation in the unknowns of Cj, of the form 

^ ' (2z(5i Pg.+ Clo.jj ~ iQo,jiPi,j2^0,i3) 

h+32+h=3 

~ 'y (2ji ~ J2 ~ 2j3)i?,i_ij4Po,j2-^o,i3 = Kij, (63) 

ii+i2+t3=t+i 

where Kij is a polynomial in the coefficients of Us<iCs, with coefficients in /3, w]. It would be 
convenient if the eight equations Eq^ j, for j G |0, 7], could define the eight unknown coefficients 
of Ci, but this is not exactly what happens, for two reasons. 

First, the equation Eq^ q involves none of the coefficients of Ci. Indeed, we see on (63), 
specialized at j — 0, that the only coefficients of Ci that Eq^ g may involve are Qi,o, Pi,o, Ri-i,o 
and But they do not occur, because Pq.o = ^ 0,1 = Qo,o = 0 (this follows from the initial 

conditions (46)). Hence this equation reads Ki^ = 0 and only involves coefficients of Us<iC 3 . 
We leave it to the reader to check that it is linear in the coefficients of Ci-i, provided i > 2. 
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Then, a similar problem happens with the sum of the eight equations Eq^ ^ for j G |0, 7]: it 
does not involve any of the coefficients of Ci either. Indeed, it reads 

y ' ~ ^Qo,jiPi,j2^0,j3) 

jl ,32 ,33 

7 

~ y (2jl ~ j2 ~ = y Ki,j- 

h,32,33 t=0 

But the left-hand side is the coefficient of in 

2Q;(<, 1)P(0, 1)7^(0, 1) - Q(0, l)p;(t, 1)Z7(0,1) - 1)P(0,1)P(0,1) 

+ R{t, 1)P'(0,1)P(0,1) + 2R{t, 1)P(0,1)P;(0,1), 

and this series is zero because P(0,1) = Q(0,1) = P^{0, 1) = 0 (see the initial conditions (46)). 
Hence this sum of equations only involves coefficients of Us<iCs. We leave it to the reader to 
check that again, this sum is linear in the coefficients of Ci_i, provided i > 2. 

These observations lead us to consider the following systeirr of eight equations: 

Xl®q*+ij,Eqz+i.0:Eqi2:Eq,_3,...,Eq,_7 \ . (64) 

3=0 ) 

Solving the system 5i gives 

Cl = (—4,8 — 2wq, 4w{q — ^) — 2/3 — 4,2/3 — 2wq] wq + 2P + 4,4rt;/3 — /3 + wq — 4,-2,wq — P — 4). 

(65) 

Moreover for i > 1, 5^ is a system of eight linear equations for the eight unknowns of Ci in terms 
of the rational functions in Us<iCs. The determinant of this linear system is 

256i^q^(3’^w {q — 4) [qu + /3^)^ , 

which is non-zero when q is an indeterminate but also when g 7 ^ 0,4 is of the form 2-1-2 cos(27r/m). 
By induction, this proves that the coefficients Pij,Qij, Rij are uniquely determined by the 
differential system and its initial conditions. Moreover these coefficients lie in <Q{q,/3,w), and 
their denominators are products of terms q, /3, w, {q — 4) and {qv + /3^). 

4.7. About possible singularities 

It remains to prove that the coefficients of the series Pj , Qj and Rj are polynomials in g, v 
and w. 

We will use three identities that we will establish later. Their proofs do not assume anything 
on the singularities of the coefficients of Pj , Qj and Rj. The first one is the characterization (49) 
of Ml, established in Section 4.8: 

2 (/32 +qv)M[ + {l + v-w (2/3 -h q)) P 3/2 -Ri = 2 + 2(3w. ( 66 ) 

The other two are established in Section 6.2: 

13 {wq -k /3) (g - 4) Qo + 9 (/3 + 2) Po + 2 {P (g - 4) {wq -k /3) t -k g) Pi = 

2/3 (g - 4) {wq - 2) {wq -k /3) t -k 2g {wq - 2), (67) 

(3 {wq -k /3) (g — 4) Qi — 2 (/3^ -k g/3 -k g) Pq — g (,3 -k 2) Pi = 

2P (g — 4) {2Pw -k wg — /3 — 2) {wq + j3)t — 2q {j3qw — 2j3w + wq — P — 2). ( 68 ) 

Let us now prove by induction on i that the coefficients of Ci have no singularity at g = 4 or g = 
—0^lv. This holds for i = 0 and i = 1. From ( 66 ), we derive that this holds for Pi_iq if it holds 
for Pi_i^ 3 , which we assume by the induction hypothesis. Now if we remove the last equation from 
the system Si (given by (64)), we obtain seven polynomial equations between the coefficients of 
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Us<iCs. Once the values of Ci are known, they are linear in Pi^, Pi^i, Pi, 3 j Qi, 0 i Qi, 0 j Ri-i.o^ 
with determinant: 

- 128z®9^/3^u; + /3q^{P'^ + 21/3^ + 36^^ + 11/3 - 6) 

+ 34/3^ + 42/3 + 12) + + 24/3 + 16) + /3^(/3 + 6)) . 

The last factor is irreducible, and this determinant contains no factor (g —4) nor ((/u + /3^). This 
proves that the coefficients in Ci are not singular at g = 4 nor q = 

Let be the equations obtained by extracting the coefficient of C in the 

equations (66), (67) and (68). The system 

7 

5] Eq,+i_^., Eq^^,,..., Eq^^,, Eq^, Eqf \ Eqf ^ 
i=o 

relates polynomially the coefficients Psj,Qs,j,Rs,j for s < z and the coefficients (in t) of the 
series For z > 1, it is linear in Pifi, ■. ■, Pi, 3 ,QifijQi,ij Ri.Oj Ri.i, with determinant 

-16^®/3®(g-4)^^;(gu + /32)^ 
and this excludes singularities at q = 0. 

Observe that the sum of (67) and (68) is divisible by /3. If we consider now the system 
{Eq,+i,o, Eq 4 ,„ ..., Eq^,,, Eq«, Eqf ^, /3-i(Eqf ^ + Eqf ^)} , 
in the same unknowns as before, we obtain the determinant 

I6i^q‘^{q - 4:)w{qiy + /3^)®, 

proving this time that the coefficients are not singular at /3 = 0. 

Finally, to rule out poles at w = 0, we resort to a different argument. First, we return to 
the case q = qm ■=‘^ + 2 cos(27r/77z) studied from Section 4.1 to Section 4.5. A first observation 
is that the series Cr involved in the invariant equation (25) are series in t with coefficients in 
Q[( 7 , u, w] . This follows from the proof of Lemma 16 in [2], using the fact that the objects denoted 
by Lr{t;y) belong to Q[q, 1/y]- A second observation is that the series Yi of Lemma 3 

satisfy Ei = 1 + f(l + j3v) + 0 {t^) where v € K*, and that their coefficients belong to K(u)[z(;]. 
This follows easily from the proof of Lemma 3, since Pj(zo) does not depend on w. Consequently, 
li := I(t,Yi) is a series in t, with a non-zero constant term that does not depend on w, and its 
other coefficients lie in K(u)[zz;]. Using these two observations, (44) now implies that P{t,x) is a 
polynomial in x with coefficients in IR(u)['u;][[f]]. By (53), the same holds for P{t,x). Similarly, 
the polynomials Q{t,x) and R{t,x) have coefficients in K(u)[zp][[t]]. 

In the next subsection, we prove (without assuming polynomiality of the coefficients of the 
series Pj, Qj and Rj), that these series, once specialized at q = q^, are indeed the coefficients 
of the polynomials P{t,x),Q{t,x) and R{t,x) constructed for q = q^. Hence, if one of these 
series had a pole at zc = 0, this would remain the case for infinitely many values q^, and this 
contradicts the fact that P{t,x),Q{t,x) and R{t,x) have coefficients in K(u)[zz;][[f]]. 

4.8. Conclusion of the proof 

We can now conclude the proof of Theorem 1. Let q be an indeterminate. We have proved in 
Section 4.6 that the differential system (47) and the initial conditions (46) define the series Pj, 
Qj and Rj uniquely as formal power series in t, and that their coefficients lie in Q[( 7 ,/3, zc]. 

Let us temporarily denote these series by Pj, Qj and Rj, to avoid confusion with the series 
denoted Pj, Qj and Rj above, which depend on a specific value of q of the form q^ := 2 + 
2 cos{2'k/ m) . Specializing the indeterminate q to qm in the series Pj, Qj and Rj gives series in t 
satisfying the differential system and its initial conditions. But we have also proved that this 
system has a unique power series solution when q = q-m- Thus Pj evaluated at q = qm coincides 
with Pj, and similar statements relate the series Qj and Rj and their barred versions. We have 
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proved in Subsection 4.5 that when q — the Potts generating function Mi is related to the 
Pj’s, Qj’s and i?j’s by (48) and (49). This means that when q = q^, 

12 (^^ + qv^ Ml + P 3 /4 + 2t (1 + — w(2f] + q)) — P 2 + 2Qo = 4t (1 + w[3j3 q)) 

and 

2 (/32 +q,y)M[ + (i + ^-u, (2/3 + q)) P 3/2 - i?i = 2 + 2/3w. 

Since all series involved in these identities have polynomial coefficients in q, and coincide for 
infinitely many values of q, they must hold for q an indeterminate. 

Let us finally prove that Mi is differentially algebraic. This follows from the uniqueness of 
the solution of our differential system in terms of power series, via an approximation theorem 
due to Denef and Lipshitz [20, Thm. 2.1]. This theorem generalizes to differential systems one of 
Artin’s approximation theorems for algebraic systems, and implies, in our context, that each of 
the series Pj, Qj and Rj is differentially algebraic. The expression of Mi = t^Mi given by (48), 
and the fact that differentially algebraic series form a ring, implies that Mi is also differentially 
algebraic. 


5. Differential system for coloured triangulations 

We now consider triangulations, and more generally near-triangulations, which are planar 
maps in which every non-root face has degree 3. We weight these maps by the number of 
vertices (variable w), the degree of the root-face (y), and by their Potts polynomial (divided by 
q). We denote by T{q, v, w, y) = T{y) the associated generating function: 

T{q, ,,,w,y) = -Y, Pm(9, , 

where the sum runs over all near-triangulations. We ignore the number of edges, which would 
be redundant: a near-triangulation with v vertices and outer degree d has 3u — d — 3 edges. In 
our first paper on coloured maps [2], we counted edges (with a variable t) rather than vertices, 
using a generating function Q{q, v, t; x, y) involving two catalytic variables x and y. ft is related 
to T{y) by: 

T{q, V, w; y) = T{y)=w Q{q, v, 0, w'^^'^y). 

Our objective is to establish a differential system for the Potts generating function of near¬ 
triangulations of outer degree 1, denoted by Ti = Ti{w). Note that this is the coefficient of y in 
T{y). More generally, we write 

Tiq, V, w,y) Td{w)y‘^, 

d>0 

hoping that no confusion arises with the mth Chebyshev polynomial T^- The root-edge of near¬ 
triangulations counted by Ti is a loop. Its deletion gives a near-triangulation of outer degree 2, 
and thus 

Ti = FT2, 

an identity that will be useful later. The expansion of Ti at order 3 reads 

Ti = i/(g -1-1- iy)w'^ -I- 1 / ((g — 1)((3' — 2 2u) i/^(q — 1 -|- j/^) 

P2v{q - 1 -h v){q -l + v'^)P u^[q - 1 -f + 0(w‘‘), 

as illustrated in Figure 3. 

Our differential system for near-triangulations is very similar to the one obtained for general 
maps (Theorem 1), but a bit simpler. 

Theorem 6. Let q be an indeterminate, P = u — I and 

D{w, x) = qv^x^ -f q)x qPv{q — 'L)w -\- P'^. 
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Figure 3. The rooted near-triangulations of outer degree 1 with v = 2 and 
V = 3 vertices, and their Potts polynomials (divided by q). 


There exists a unique triple {P(w,x),Q{w,x), R{w,x)) of polynomials in x with coefficients in 
Q[( 3 ', u] [[w]], having degree 3,2 and 1 respectively in x, such that 


and 


x^]P{w, x) 

= 1 , 

P{G,x) = 
Q( 0 ,x) = 

x‘^{x 1/4), 

x{2h'x 1 ), 

(69) 

1 d , 

f R^ \ 

1 d / 


(70) 

Rdx ' 

yPD'^ j 

Q dw \ 

PD^) ■ 


Let Pj(w) = Pj (resp. Qj, Rj) denote the coefficient ofx^ inP(w,x) (resp. Q(w,x), R{w,x)). 
The Potts generating function of near-triangulations of outer degree 1, denoted by Ti{w) = Ti, 
can be expressed in terms of these series using 


20 n^qTi - -h 4z/Qo + {Qi - 1 ) (Qi + i' - 3) + 2 v {qv - 24^ -&q)w = 0 . (71) 


An alternative characterization of Ti is 

2izqT[ =Ri- q{P - 1) -f 8/3. 


(72) 


The series Ti is differentially algebraic, that is, satisfies a non-trivial differential equation with 
respect to the vertex variable w. The same holds for each series Pj, Qj and Rj. 


Eq. (70) is, in compact form, a system of eight differential equations in w relating the 9 series 
Pj, Qj, Rj. Since P 3 is given explicitly, we have in fact as many equations as unknown series. We 
will see that Q 2 = 2u. In fact, we will prove that (70), combined with the initial conditions (69), 
determines uniquely all series Pj, Qj, and Rj. For instance, 

Pq = -Pw -T P ( 8 q-i- q{q - 12)/3/4 - {q-\- 6)/3^ - 3/3^) w'^ + 0 {w^). 

We expect Ti to satisfy a differential equation of order 4 (Section 6.2). We will work out in 
details several special cases in which Ti satisfies a second order DE (Sections 7 to 9) . 

The proof of the above theorem is similar to that of Theorem 1, and we mostly tell where 
these proofs differ, without giving details otherwise. 


5.1. An equation with one catalytic variable 

As in the case of general planar maps, we begin with an equation for the series T{y), taken 
from [2]. We assume that <7 = 2-1- 2cos(2/c7r/m) with k and m coprime and 0 < k < 2m. We 
still write /3 = u — 1. We introduce the following notation: 

• I{w, y) = I{q, V, w; y) is a variant of the generating function T{q, v, w; y): 

I{w,y) = yqT{q,v,w,y) - ^ 

y y^ 

• N{y,x) and D{w,x) are the following (Laurent) polynomials: 

N{y, x) = P(4 - q)y -k qvx + P{q - 2), 

with y = 1 /y, and 

D{w, x) = qv^x^ /3(4/3 -f q)x qPv{q — 4)w -\- P'^. 

We still denote by Tm the mth Chebyshev polynomial of the first kind. 


( 73 ) 
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Proposition 7. There exist m + 1 formal power series in w with coefficients in Q{q, v), denoted 
Cq(w), ... ,Cm{w), such that 

D{w, I{w, y)r/^ = f2 Cr{w)nw, yy. (74) 

\2^yD{w,I{w,y)) J ^ 

Proof. This follows from Corollary 12 and Lemma 19 of [2]. One must prove that the series in 
t denoted Cr in [2] (which coincide with the series Cr of the present paper, with t = are 

not only series in t but in fact series in y. This can be done by following carefully the proof of 
Lemma 19 of [2], using the fact that the series denoted K{ty) is a series in t^. This implies that 
the series denoted PSij and PTij are also series in t^, and one concludes using Eq. (93) of [2]. 


By expanding the invariant equation (74) near y = 0, we can express the series Cr in terms of 
the derivatives of T with respect to y, evaluated at y = 0. We will need the following expressions 
of Cm, ■ ■ ■, Cm- 3 - The first three are explicit in terms of q, v and w, but the last one involves 
the series Ty{0) = Ti. 

Lemma 8. Denote jd = v — 1, and recall that q = 2 + 2 cos{2kTT/m). We have: 

{-y^c„ 


(gjy2)m/2 


= 1 , 


(-1)'=C„1 

m/3 

{^qyl'ynl’l-l 

2 

{-\yCm-2 

m/3 

(^qjy2'^m/2-2 

24 

( 

{-If Cm-3 

1 


(4/3 + q + I3m{q - 4)), 

^ = -yy ( 12q^iy^w{q — 4) + 48^^/3(m — 1) + /3(m — 1)((7 — 4)x 

(dq + 12 + Qj3{mq + 4) + I3‘^{imfq — Amf + mq + 20m — 12)^ ^ , 


,2„3,,4o2 , 


(gj,2)m/2-3 - 2 ^ ^ 

m{m — l)/3^ 
720 


180q^v^w{q — 4) (^{q — A)j3m + 8/3 + 2y^ + 960^1/^ (m — 2) 


+ 720/3^^(m — 2){q — 4)(/3m — /3 + 2) + jd{rn — 2){q — 4)^ (^0^mf{q — 4) 

+ + 16/3 + 5y) + Pm{2(d^q - 8/3^ + 15/3y + 360/3 + 60q + 180) 

- 180/3^ + 360/3 + 60g + 300^ 

where Ti is the Potts generating function for near-triangulations with outer degree 1. 

Remark. When comparing this lemma with its counterpart for general planar maps, Lemma 2, 
we observe that we have no term Tm{l) nor T(,^( 7 ) here. This is because 

N{y,I{w,y)) ^/q n / \ 

7 = lim- , = = cosi/cTT/m), 

y^0 2^D{w,Iiw,y)) 2 

so that Tmil) = (—1)^ and T(„( 7 ) = 0. 


Proof. We multiply the invariant equation (74) by so that it becomes not singular at y = 0, 
and expand it around y = 0. We extract successively the coefficients of y^,y‘^,y‘^ and y® and 
obtain in this way expressions of Cm, ■ ■ ■, Cm -3 (the coefficients of odd powers of y do not give 
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more information). For instance, since 

y^I{w,y) = l + 0{y), 
y'^N{y,I{w,y)) =qiy + 0{y), 
y‘^D{wJ{w,y)) = qv"^ +0{y), 
extracting the coefficient of gives 

and the expression of Cm follows since Tm(Y^/2) = (—1)^- As we extract the coefficients of 
higher powers of y, derivatives of T^, taken at the point y^/2 = cos{kn/m), occur. We can find 
explicit expressions for them in terms of q, k and m using T),^(^/g/ 2 ) = 0 and the differential 
equation 

(l-u2)T"(u)-uT'^(u)+m2T„(u) =0. 

Similarly, derivatives of y'^I(w, y) with respect to y occur, and this is why the expression of Cm-s 
involves the series Ti. B 


5.2. Some special values of y 


We now work out the counterpart of Lemma 3. We still denote n = [m/2j. 


Lemma 9. Let q = 2 + 2 cos(27r/TO), and denote fd = v — 1. There exist m — 2 distinct formal 
power series in w, denoted Yi,..., that have constant term 1 + 0(/3) as /3 —>■ 0 and satisfy 

N{Y,I{w,Y)Y = 'icos{jTT/m)^D(w,I{w,Y)), (75) 

for some j S |l,u]. Their coefficients are algebraic functions of (3 ouer K. 

Let us denote Ii{w) := L(w,Yi). This is a formal power series in w with coefficients in 
the algebraic closure o/IR(/3), denoted by M(/3). The (m — 2) series L are distinct, and for 
1 < i < m — 2, 

dl 

L(0)^{0,-l/4,-l/(2u)}, —{w,Yi)^Q. 

The key difference with Lemma 3 is that the coefficients of the series Yi do not have their 
coefficients in M(/3), but in its algebraic closure. Let us illustrate this by an example. 

Example: bicoloured triangulations. We take q = 2, that is, to = 4 (and n = 2). The 
relevant values of j are 1 and 2, and we want to solve 

N{Y,I{w,Y)f = 2D{w,I{w,Y)) 


for j = 1, and 


N{Y,I{w,Y)) = 0 

for j = 2. Using the definitions of N and D, this reads, for j = 1, 
4y3 (2/3r -2(3+ Y -2) T{Y) - 4Y^ (1 + /3) w + (F - 2) - 


2^F + 2/3 - 2F + 2) = 0, 


and for j = 2, 


2Y^ (1 + /3) r(F) + /3 + 1 - F = 0. 


Recall that we are interested in series Y{w) with constant term 1 + 0{(3). 
solution for each of the above equation, which satisfy 


We find one such 


Fi(i(;) 

Y2{w) 


/3 + 1- 


/3 

1 + (3 3- 0(w). 


0(w), 


□ 


Proof. We denote qj = Acos{jTr/m)^ = 2 + 2cos{2jTr/m). Observe that q = qi. We assume that 
TO = 2n + 1 is odd. The proof is similar in the even case. 
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Existence of the series 1^. We are looking for series V = Y(w) solutions of (75). Multiply¬ 
ing (75) by we obtain the equation $j(w,V) = 0, where 

$j (u;,y) = y‘^N{yJ{w,y) f - qjy^D{w, I{w,y)). 

Observe that y^I{w, y) is a formal power series in w with coefficients in y]. Hence the same 
holds for y‘^N{y,I{w,y)) and y'^D{w, I{w,y)), and finally for ^j{w,y). Expanding I, N and D 
at first order in w gives: 

y^I{w,y) = l-y + 0{w), (76) 

y^N{y,I{w, y)) = qv - y{q + 2/3(9 - 2)) + (g - 2)/3y^ -h 0{w), 

y^D{w,I{w, y)) = qv"^ - 2 qv'^y + {q + 39/3 + [q + 4)/3^) ?/^ - /3 {A/3 + q)y^ + /3‘^y'^ + 0{w). 

(77) 

The equation A>j{w,Y) = 0 thus reads 

Pj{Y) + wSj{w,Y) =0, (78) 

where 

Pjiy) = {yv - y{<i + 2^(9 - 2 )) + (9 - 

- Q] - 2 qy'^y + (g + 39/3 -k (9 + A)0^) y"^ - p (4/3 + q)y^ + /3^/) 

and Sj{w,y) is a power series in w with coefficients in y]. In particular, the coefficient of 

in Y, denoted yo, must satisfy Pj{yo) = 0. The roots of this quartic polynomial can be seen 
as Puiseux series in P with coefficients in C (see for instance [39, Ch. 6]). Let us focus on the 
roots that are finite at /3 = 0 and have constant term 1. Using Newton’s polygon method, we 
find that they read j/o = 1 + /3u -I- 0(/3^), where v must satisfy the equation (34) that we studied 
when constructing the series Yi for general planar maps. Thus when j > 1 we find for Pj two 
distinct roots with constant term 1, of the form y^ = 1 + Pv^ + 0{P^), and only one such root 

2/0 = 1 + Pvi + 0{P^) when j = 1 (with = l/q). 

Now having fixed one root yo of Pj, let us return to (78) and extract from it the coefficient 
of w^, for p > 1. This gives 

Pj{yo)[w^]Y -f (expression involving only [w^]Y for i < p) =Q. 

Since yo is not a double root of Pj, this allows us to compute the coefficient of in Y by 
induction on p, and thus determines Y completely. Moreover, is an algebraic function of 

P over M. 

As argued in the proof of Lemma 3, the (m — 2) values ,... ,v^,V 2 , ■ ■ ■ ,v~ are distinct. 
They give rise to (m — 2) distinct series Y{w) satisfying +(0) = 1 + 0{P) and (75) for some 
j S |1, n]. We denote them Yi,..., y „_2 with ^(0) = 1 + Pvi + 0{P^), where as before, 

{wi,...,U„_ 2 } = {vt,...,V+,V 2 ,...,V-}. 

Clearly these series are non-zero, and we have proved the first statement of the lemma. 

Properties of Yi and Jj := I{w, Yi). Let us prove that the series A are distinct. Returning 
to (76) gives 

/i(0) = -Pvi + 0(/3^). 

This proves that the series A are distinct (since the vPs are distinct), and also that li{0) is 
distinct from 0,—1/4 and —l/{2v). 

Let us now prove that D{w, Ii{w)) ^ 0. Returning to (77), we find that: 

D(0,li(0)) = P^(qv^ - qvi + 1) + 0(/3^). 

Comparing with (35) shows that this is non-zero. 

We finally check that 


— (0,y(0)) = -l + 4/3u, + 0(/32) 
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is non-zero, as stated in the lemma. a 

5.3. Some polynomials with common roots 

We still assume that q — 2 + 2cos(27r/m). We denote as before: 

ra 

C{w, X) = ^ Cr(w) X''. 

r—0 

Thus C{w,x) is a polynomial in x with coefficients in IR(u)[['u;]]. This is also true of D{w,x) 
(see (73)). The coefficients of D are explicit, but those of C are unknown. Here is now the 
counterpart of Lemma 4. 

Lemma 10. Each of the series li := I{w, Yi) defined hy Lemma 9 satisfies 

= D{w,Ur, 

m f)r) 

D{w,L) — {w,Ii) = — 

m f)r) 

D{w,L) — {w,Ii) = -C{w,h) — {w,h). 

The first two identities imply that L is actually a double root of . 

The proof is identical to that of Lemma 4, with the variable t replaced by w. 

Proposition 11. Let li := I{w,Yi) be the (m — 2) series defined in Lemma 9. There exists a 
triple {P{w, x), Q(w, x), R{w, x)) of polynomials in x with coefficients in M(/3)[[ry]], having degree 
at most 3, 2 and 1 respectively in x, such that 

m—2 

C{w, x)^ — D{w, a;)"* = P{w, x) (a: — , 

i=l 

m-2 

D{w, x)C'^{w, x) - —D'^{w, x)C{w, x) = Q{w, x) (a: - R), 

i=l 

m-2 

D{w, x)C'^{w, x) - —D'^(w, x)C(w, x) = R{w, x) (a; - R). 

i=l 

Proof. The proof is almost the same as that of Proposition 5, with the variable t replaced 
by w. The only difference is that here the polynomial C{w,x)'^ — D{w,x)'^ has degree at most 
2m — 1 (instead of 2m), and the polynomial D{w,x)C!^{w,x) — 'fD'^{w,x)C{w,x) has degree 
at most m — 1 (instead of m). This is because the coefficient of a:^™ in C{w,x)^ and D{w,x)'^ 
is (by Lemma 8), and the coefficient of x"^ in D{w,x)C!^{w,x) and D'.^{w, x)C{w, x) 

is — (again, by Lemma 8). B 

5.4. Differential system 

We still assume that q = 2 + 2cos(27r/m). 

A differential equation relating P, Q and R. Starting from Proposition 11, one first proves 
that 

Qdw ypD^) ~ Rdx ypD^) ■ 

This argument is the same as in Section 4.4. The above equation coincides with (70), but 
with hats over the letters P, Q and R. The series P, Q and R occurring in Theorem 6 will be 
normalizations of P, Q and R by multiplicative constants, independent from w and x. 
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The leading coefficients of P,Q,R. Proposition 11 gives 
[x^]P{w,x) = {Ciw,xy - Diw,x)^) 

[x'^]Q{w,x) = [x™] (^D{w,x)C'^{w,x) 

[x^]i?(w,x) = [x™”^] (^D{w,x)C'^{w,x) 


- -^D'^{'w,x)C{w,x)j , 


m _ 


- -^D'^iw,x)C{w,x)j . 

Recall that D is given explicitly by (73), and that Lemma 8 gives the leading coefficients of 
C{w,x). This allows to us determine the leading coefficients of P,Q, R: 


P3 := [x^]P(ri;, x) = I 3 ‘^{q — 4)g™ ij/2m 2 

Q 2 := [x^]Q{w,x) = 

Xal 


f? r ™ / 3 % m /2 m T,. , 9(/3-1)-8/3^ 

Pi := [x ]P(w, x) = —^—(g - 4)g ' v \qT^ + - - ^ -- I . 


(79) 


Let us define 


P{w,x) 2i^Q{w,x) 2i^R{w,x) 

P{w,x) = —-, Q{w,x) = -- and R{w^x) = 


(80) 

(81) 


P 3 Q 2 Q 2 

Then P, Q, R satisfy (70) and moreover, 

[x^]P(ic, x) = 1, [x^]Q{w,x) = 2u, [x^]R{w,x) = 2vqT[+q{l3 — 1) — ij3. 

In particular, the first of the initial conditions (69) holds. (We do not give Q 2 explicitly in 
the statement of the theorem, as its value follows from the differential system and the initial 
conditions.) 

The case tn = 0. It remains to determine the values P(0,x) and <3(0, x). Proposition 11 gives 


m 

-2^ 


m—2 

p(o,x) 

2^1 

(82) 

m-2 

Q(0,x) (x- /i(0)). 

(83) 


i=l 


We will now combine our knowledge of C{Q, x) obtained in [2] with the properties of the series 
Ii(w) gathered in Lemma 9 to determine P(0,x) and Q(0,x). 

The series (7(0,x) can be determined using results in [2, Lemma 19] (with t = There, 

it is proved that 

C( 0 ,x 2 - x) ^ ( 27 ^) ’ 

J\f = iV(l/x, x^ — x) = /3(4 — q)x + qi 2 {x'^ — x) + I3{q — 2) 


where 

and 


V = P(0, x^ — x) = q{l + /3)^(x^ — x)^ + /3(4/3 + q){x^ — x) + /3^. 
Recall the expression (58) of — 1. Thus (??) gives 

L^J 


1 


(7(0,x2-x)2-P(0,x2-x)™ =-(Af2-4P)(Af^) 


i=i 


'cos 


2 


(84) 


m—2 


= P(0, x^ - x) (x^ - X - /i(0))" 


i=l 
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by (82). We want to prove that x^{x + 1/4) is a factor of P{0,x), or equivalently, that x‘^{x — 
^Y{x — 1/2)^ is a factor of P(0,a;^ — x). Since by Lemma 9, li(0) ^ {0, —1/4}, it suffices to 
prove that x‘^{x— l)'^{x— 1/2)^ divides the right-hand side of (84). Using the above values of Af 
and V, we first find a factor (a; — 1)^ in {Af^ — 41?). Then, using 4cos^(7r/m) = q, we find that 
the term obtained for j = 1 contains a factor {x— 1/2)^. Finally, using 4cos^(27r/m) = {q — 2)^, 
we find that the term obtained for j = 2 contains a factor a;^, provided m > 5. When m = 3, 
that is, (J = 1, a factor x^ is found in the term obtained for j = 1. When m = 4, that is g = 2, a 
factor x^ comes out from the term . 

We now wish to prove that x{2vx + 1) divides Q{0,x). We have just seen that a; = 0 is a 
double root of P{0,x). Using the same argument as in Section 4.4, this implies that x divides 
Q{0,x). For the remaining factor, we have to prove that (2v{x^ — a:) -I- 1) divides Q{0,x'^ — x). 
Let us combine (83) (taken at x'^ — x) with (??). We obtain: 

m—2 

(5(0, x^ — x) (a;^ — X — li(0)) = P(0, x^ — x)C^(0, x^ — x) — ^Z?(,(0, x^ — x)C(0, x^ — x) 

i=l 

■pl+m /2 g / J\f \ / Af \ 

~2x-ldx[2VvJ "‘\2VvJ' 

A factor (2u(a;^ — a;) -I- 1) is found in ^ Lemma 9 ensures that /i(0) ^ —l/(2u), and 

thus this must be a factor of Q. 

5.5. The Potts generating function of triangulations 

It follows from the first two identities of Proposition 11 that 

P (dc, - jKcy = . 

Using the normalisation (80), this gives 

P (pc' - ^D'^cy = ?^(g/4- l)g2 (C^ - P-), (85) 

since Q 2 = qm?j3'^{l + ^)^(g/4 — 1 )P 3 by (79). 

Recall that the leading coefficients of C are given in Lemma 8 , and that C^ — P"* has degree 
2m — 1. Recall also the known values (81) of P 3 and (52- Extracting the coefficients of 
^ 2 m+ 2 ^ and in (85) first gives (for the coefficient of x^’”’*'^) a tautology, then (for the 

coefficient of x^'"^^) an interesting relation between P 2 and Qi, namely 

4uP2 = 4Qi -t u - 4, (86) 

and finally (for the coefficient of x^’”’*'^) and expression of Pi which we transform into (71) 
thanks to ( 86 ). 

The expression (72) of was obtained in (81). 

5.6. Uniqueness of the solution 

We now want to show that the differential system of Theorem 6 , together with its initial 
conditions, uniquely defines the nine series Pq, Pi, P 2 , P 3 , (5oi Qh Q 2 , R 01 Pii whether q is 
an indeterminate or a real number of the form 2 -|- 2cos(27r/m). The proof parallels the case of 
general planar maps given in Section 4.6 (with t replaced by w), and is actually a bit simpler. 
We merely sketch the main steps. 

The system can be written as in (62): 

2Q'^PD - QP'^D - 2QPP/, = 2R!^PD - RP'^D - 2PPP(,. (87) 

Both sides are polynomials of degree at most 7 in x. Recall that P 3 = 1. Extracting the coefficient 
of x^ gives Q 2 {w) = 0, which, with the initial condition ( 52 ( 0 ) = 2u, implies Q 2 {w) = 2v. We 
are left with seven equations and seven unknown series. 
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We denote Pij = [w^]Pj, and similarly for Q and R. Let Ci = {Pi^, Pt^i, Pi^., Qi,o,QiX, Ri-i, 0 t 
Ri-i^i). The right-hand side of (69) gives 

Co = (0,0,1/4; 0,1; 0,0). 

The next steps are the same as in Section 4.6. For i > 1 and j G |0, 6 ], we denote by Eq, j the 
equation obtained by extracting the coefficient of in (87). Again, the equation Eq^ g only 

involves series from Us<iCi because Pq.o = ^o,i = Qo,o as in the case of general maps. However, 
we do not have here the second difficulty that came from the factor {x — 1) in P{0,x), P^{0,x) 
and ( 5 ( 0 , a:), and we consider the system 

Si = {Eq,+io,Eq,_i,Eqi_ 2 ,---,Eq,_g} . 

Solving the system for the unknowns of C\ gives 

Cl = (-^, ^2 + ;3g/2-4/3, 4 / 32 + 2/3g-2(7; ^0^X2qf+ -2q-, -2/3, Pq-S^-q). 

Then for z > 1, 5^ is a system of seven linear equations for the seven unknowns in Ci, the 
determinant of which is: 

l^q^0 {q - 4) (/3 + 1)" (/3 - 1)" ( 4/32 - q) /2. 

We conclude as in Section 4.6. 

From this point on, we argue as in Section 4.8 to prove that, for q an indeterminate, the series 
Pj, Qj and Rj are differentially algebraic and satisfy (71) and (72). It then follows from (71) 
that Ti is differentially algebraic too. 

5.7. About possible singularities 

It remains to prove that the coefficients of the series Pj, Qj, Rj are polynomials in q and v. 
First, we note that (72) implies that this is true for i?i. Consider the system of (only) six 
equations: 

{Eq^-i-i Q, Eq^ 2; ■ • ■; Eq^ g}. 

It relates polynomially the coefficients in Us<iCs. For z > 1, it is linear in Pi^, Pi^i,Pi^ 2 iQi,o^ 
Qi^i, Ri-ifi once the values Pi^ and Rq^ are known, with determinant: 

-l^q00 - 1) ((?" - ^q00 -P + l) + 2q0 - lf{0 - 4/3 + 1) - - 1 )^) . 

This determinant contains no factor (g —4), (/3-I-1) nor {A0 — q). By induction on i, we conclude 
that the denominators of the coefficients of Pj, Qj and Rj consist of factors g,/3 and (/3 — 1). 

To show that these factors do not occur, we use two identites that we will prove in Section 6.1 
(without assuming polynomiality of the coefficients...): 

iyq{q - 4)Qo - (4^ -k q)Ro + 2{vqw{q - 4) -k j3)Ri = 2/3{q - 4){vqw{q - 4) -f /3) 

and 

p/3g(g - 4)Qi - 2zz^gi?o -k /3(4^ -k q)Ri = 2{q - A)0{AP + q). 

Let Eq( and Eq/ be obtained by extracting the coefficient of w® in these equations. Consider 
the system 

{Eq,+i 0 , Eq, 1 , Eq,_ 2 , Eq^ 3 , Eq', Eq"}. 

It relates polynomially the coefficients Psj,Qs,j and Rg^j for s < i. For z > 1, it is linear in 
Pi,o,Pz,i,Pi, 2 ,Qi,o,Qi,i^Ri,o, with determinant 

-2zV°(/3 + l)(9-4)(4/3 + g). 

We have thus ruled out factors (/3 — 1). To rule out the factor q, we note that Eq/ is a multiple 
of q (because by (72), all coefficients for z > 0, are multiples of q). So if we replace the last 
equation of the previous system by Eq/ /q, the factor q disappears from the determinant. 

We are left with the factor j3. This time we form the system 

{Eq*, 3 , Eq,_ 4 , Eq^ g, Eq^ g, Eq', Eq/}, 
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which is linear in the unknowns Pi^, Pi^i, Pi, 2 i Qi, 0 ) Qi,i) Ri, 0 i with determinant 

-2i^q^{q - 4)(/3 - l)(/3 + 1)1°(4/32 - q). 

There is no factor /3, and we conclude that the coefficients of Pj, Qj and Rj belong to Q[( 3 ', u]. 


6. Simplifying the system: five non-differential equations 

In this section, we derive five equations between the unknown series Pj, Qj and Rj, which do 
not involve their derivatives (however, one of them involves the derivative of Mi or Ti, depending 
on the family of maps we consider). 

Since the system obtained for triangulations is a bit simpler than for general maps (the degrees 
of P and R being smaller), we begin with this case. 


6.1. Triangulations 


We start from the system of Theorem 6 . It involves nine series in t denoted Pq, ..., P^, 
Qo, ..., Q 2 , Rq, Ri- One of them is given explicitly in the theorem: P3 = 1 . Moreover, we have 
derived from the system in Section 5.6 that Q2 = 2u. This leaves us with seven unknown series, 
related by seven differential equations. These equations are the coefficients oi ,... ,x^ in (87). 

We are going to derive three linear (non-differential) relations between the unknown series by 
letting X tend to 00 or to one of the two roots (in x) of D{w,x). Extracting the coefficient of x® 
from (87) gives 

Q'l = 

Using the initial conditions (69), we obtain 

UP 2 = Qi + u/4-1. ( 88 ) 

This identity was already obtained by expanding (85) in the proof of Theorem 6 (see ( 86 )). 


Let us now specialize (87) at the two roots Si{w) and 62 {w) of D{w,x). These roots, seen as 
series in w, satisfy: 




0 (w), 


with (3 = V — 1 as above. Since P(0,x) = x^(x -I- 1/4), the Si's are not roots of P, and thus 
taking the limit x —>■ <5^ in (87) gives, for i = 1,2, 


Q{w, 5^{w))D',^{w, 5i{w)) = R{w, Si{w))D'„^{w, Si{w)). 


Since D and its roots are explicit, this is a linear system in Qq, Qi, Rq and Pi, symmetric in di 
and 62 . Solving for Qo and Qi gives 


uq{q - 4)Qo - (4/3 -k q)Ro + 2{vqw{q - 4) -k /3)Pi = 2/3(g - P){uqw{q - 4) -f /3) (89) 


and 

u/3q(g - 4)(5i - 2u^gPo-I-/3(4/3-I-g)Pi = 2(g - 4)/3^(4/3-|-q). (90) 

Let us finally recall the two characterizations of Ti obtained in Theorem 6 : 

2Qv‘^qTi - du^Pi -k 4uQo + (Qi - 1) (Qi + u - 3) -k 2v {qu - 24/3 -Qq)w = 0. (91) 


and 

2vqT[ = Pi - g(/3 - 1) -k 8/3. (92) 

We call Eqs. (88-92) the five generic non-differential equations: they involve the derivative T[, 
but no derivative of the other unknown series Pj, Qj and Rj. And of course they hold for 
indeterminates q and v. 


Prediction of the order of Ti. Using (in this order), equations ( 88 ), (92), (90), (89) and (91), 
we obtain expressions of Qi, Pi, Rq, Qo and finally Pi as linear combinations of P 2 , Pi and T[ 
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(with coefficients in Q{q,i/,w)). We only give the connection between Qq, P 2 and Ti, which we 
will use later: 

2 [Av'^q^w + - q)) T^jv + 2qi^Qo - P{q + 4 / 3 )P 2 = (4^5 + q){Aqvw - /3/4). (93) 

Recall that Ti = VT 2 where T 2 counts near-triangulations with outer degree 2. 

We have thus got rid of five unknown series, but introduced a new one, namely Ti. Reporting 
the expressions of Qi, i?i, i? 0 : Qo and Pi in our differential system gives a system of three equa¬ 
tions in Pq j P 2 and Ti, of order 1 in Pq and P 2 and order 2 in Pi. Thus the field generated over 
Q(( 7 , V, w) by Pq, P 2 , Ti and their derivatives has transcendence degree at most 4, and we expect 
Ti to satisfy a differential equation of order at most 4. We have been able to eliminate Pq and 
to obtain two differential equations in P 2 and Ti : one of order 1 in P 2 and 2 in Pi, the other of 
respective orders 2 and 3. Going further with the elimination seems to require heavy computer 
algebra, and we have failed to obtain an explicit equation for Pi. We study in Sections 7 to 9 
three special cases, of combinatorial interest, where the order of Pi is only 2 . 

6.2. General planar maps 

We now proceed similarly with the differential system of Theorem 1. ft involves eleven series 
in t denoted Pq, ..., P 4 , Qq, ..., Q 2 , Rq, ..., R 2 - Two of them are given explicitly: P 4 = 1 and 
R 2 = 1 ^ + 1 — w{q + 2/3). We have also seen in Section 4.6 that Q 2 = 1. This leaves us with eight 
unknown series, related by a system of eight differential equations. 

As in the case of triangulations, we can derive three linear (non-differential) relations between 
them by considering the system (62) as x approaches 00 , or one of the two roots (in x) of the 
quadratic polynomial D{t,x). The derivation is the same as in the previous subsection, and we 
simply give the resulting three identities. 

The first one was already obtained in the proof of Theorem 1 by expanding (60) in x (see (61)): 

P 3 = 2Qi+At{l + iy)-4tw{2(3 + q). (94) 

The other two are the counterparts of (89) and (90), and they read: 

13 {wq + P) {q - 4) Qo + q (13 + 2) Rq + 2 {(3 {q - 4) {wq + I3)t + q) Ri ^ 

2/3 {q - 4) {wq - 2) {wq + I3)t + 2q {wq - 2), (95) 

13 {wq + I3){q- 4)Qi-2 {j3^ + q/3 + q) Rq - q {f3 + 2) Ri = 

2(3 {q — 4) {2l3w + wq — 13 — 2) {wq + l3)t — 2q {/3qw — 2l3w + wq — (3 — 2). (96) 

Let us finally recall the two characterizations of Mi = t^Mi obtained in Theorem 1: 

12 (/3^ -f (/p) Ml -f P 3 /4-|-2t(l-|-j^ — w{2j3 q)) P 3 — P 2 -|- 2Qq = 4t (1 - 1 - ui(3/3 -|- q )), (97) 

and 

2 (/32 + q,y)M[ + {l + , 2 -w (2/3 -k q)) P 3/2 - Pi = 2 -k 2/3w. (98) 

Again, we call Eqs. (94-98) the five generic non-differential equations: they involve the derivative 
M'l , but no derivative of the other unknown series Pj , Qj and Rj. And of course they hold when 
q, V and w are indeterminates. 

Prediction of the order of 7\Pi. Using (in this order) the five equations (94), (98), (96), (95), 
and (97), we can express Qi, Pi, Rq, Qq, and finally P 2 as polynomials in P 3 , Mi and (with 
coefficients in Q{q,iy,w,t)). We thus get rid of five unknown series (but introduce a new one, 
namely Mi). Reporting these five expressions in our differential system gives a system of four 
equations in Pq,Pi,Pq and Mi, of order 1 in Pq,Pi,Pq and order 2 in Mi. We thus expect Mi 
to satisfy a differential equation of order at most 5. We study in Sections 8 and 10 two special 
cases, of combinatorial or physical interest, where the order is only 3. 
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7. Properly coloured triangulations: u = 0 


We specialize the differential system obtained for triangnlations to u = 0, and recover the 
solution of the problem studied by Tutte between 1973 and 1984. Recall that, for generic values 
of q and V, we have Ti = uT 2 , where Ti counts g-coloured triangnlations with outer degree i. 
Thus Ti vanishes when u = 0, and we focus on T 2 instead, as Tutte did. 

Theorem 12 (Tutte [50]). Let T 2 = T 2 {q,w) denote the generating function of properly q- 
coloured near-triangulations of outer degree 2 , counted by the number of vertices (as before, the 
root-vertex is coloured in a prescribed colour). This series is characterized by 

2(1 - q)w + iw + IOT 2 - 6wT^)T^' + (4 - q){20T2 - 18wT^ + 9w‘^T^) = 0, 

with the initial conditions T 2 = 0{w^). Equivalently, T 2 = J2n>2'^niq)w^, where a„ = an{q) is 
given by 02 = q — 1 and for n > 0, 

n 

{n + l)(n + 2 )a „+2 = {q — 4)(3n — l)(3n — 2)a„+i + 2^ i(z + l)(3n — 3z + l)ai+ia„+ 2 -i. 

i=l 

Remark. Take a properly coloured near-triangulation of outer degree 2, having at least one 
finite face (that is, having at least three vertices). If we delete its root-edge, we obtain a properly 
coloured triangulation with the same number of vertices. This means that when u = 0, 

T2 = (g - l)w'^ -\- T3, 

where T 3 counts properly g-coloured triangnlations (with the root-vertex coloured in a prescribed 
colour). The series H occurring in (1) is qT 2 , and Tutte’s differential equation (1) is equivalent 
to the above theorem. 


Proof. We specialize Theorem 6 to u = 0, that is /3 = —1. The polynomial D{w,x) becomes 
independent of w and has degree one in x'. 

D{w,x) = {A —q)x-\-\. (99) 

Let us now specialize the five generic non-differential identities (88-92). The first three give 

Qi = 1, Rq = 2 and i?i = 2(4 — q). (100) 


Observe that R(w, x) = 2D(w, x). The identities (91) and (92) then become tautologies at u = 0, 
and we still need to relate T 2 to the series Qj and Pj. We use for that the equation (93) . Once 
specialized at u = 0 , it reads: 

2r^-HP2 = l/4. (101) 

This will replace the identity (92) relating T{ and i?i. We proceed without deriving a counterpart 
of (91). 

We now go back to the differential system (70), which simplifies into 


d fl 


dx \P 


1 d 


2 ^ -] = ^ 


Q 2 

P 


( 102 ) 


QD dw 

This is a rational expression in x. Given the form (99) of D, we find convenient to expand the 
numerator of (102) in powers of X = a; -|- 1/(4 — g) rather than x. We write accordingly: 

P(w,x)=X^ + P 2 X^ + PiX + Po, Qiw,x) = X + Qo, R(w,x) = 2(4 - q)X, 

where the Pj’s and Qq are series in w with rational coefficients in g. Since we have taken 
into account the values (100) of Qi, Rq and i?i, we have only four unknown series left. The 
identity ( 101 ) gives 

P = g + ^ _ 9T' 

' 4(g-4) 

The following initial conditions are translated from (69): 

9 f, 2 -kg 


(103) 


^^o(O) = 


4(g-4)3’ 


Pi(0) = 


4(g-4)2 


^2(0) = 


8 -kg 

liq^y- 


Qo(0) = 


1 

g^- 


(104) 
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Expanding in X the numerator of (102) gives a system of four differential equations between the 
four series Po,Pi, P 2 and Qq: 


2PoQo - Qo^o = 0, 

2 Q 0 P 1 - QoP'i -Pi + 2(4 - q) A = 0, 

2Qo A - QoA - a + 4(4 - g) A = 0, 

2Q' - A + 6(4 - g) = 0. 

The first and last equations are readily solved. With the initial conditions (104), the last one 
gives 


Qo — 2 A + 3(g — 4)ui — 
= —A + 3(g — A)w + 


8(g-4)’ 

1 


g-4’ 


thanks to (103), while the first one gives 


4(g - 4)Po = qQl = q[ "A + 3(g - A)w + 


q-A 


(105) 


(106) 


Reporting the expressions (103) and (105) of A and Qo in the third differential equation gives 
an equation in Pi and T 2 that we can integrate (using the initial conditions (104)): 


A = (A)' + A 


^6(g — A)w 


g + 4 \ 

2(g-4); 


w 

+ 10(4-g)r2 + -(g + 8) + 


2 + ^ 
2(g-4)2- 


Reporting this, as well as (105) and (106) in the second differential equation finally gives Tutte’s 
equation for the generating function of properly g-coloured near-triangulations of outer degree 2. 
The rest of the proof is straightforward. B 


8. Four colours: q — A 

Several signs suggest that the case of four colours should be simpler. First, this is visible on 
the differential equation and recurrence relation of Theorem 12. Then, the polynomial D{t,x), 
given by (45) (or its counterpart D{w,x) given by (73) for triangulations) becomes a square 
when g = 4, independent of the size variable t (or w). We use this simplification to derive a 
differential equation for the generating function of four-coloured maps, of a relatively small order 
(3 for general maps, 2 for triangulations). 


8.1. Triangulations 

Theorem 13. Let Ti = Ti(i',w) be the generating function of four-coloured near-triangulations 
of outer degree 1, counted by the number of vertices (w) and the number of monochromatic edges 
(ly) (as before, the root-vertex is coloured in a prescribed colour). This series is characterized by 
the initial conditions Ti = vifi -\- v)uP' -\- 0{w^) and the following differential equation of order 2 
and degree 6, which we write in terms of S = 2Ti — w: 


P{S,S',S”) = 0 
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with 


P{X, Y, Z) = + iP^awY'^Z'^ + l02A0iy'^p'^w^Y^Z^ - 5P^aXY^Z^ 

- 8P^w{10Xl'^ - a^)Y^Z^ + 768v^aP^w'^Y^Z"^ + 58982Aiy^w'^Y^Z^ + A9152L''^Paw^Y^Z^ 

- IbP^a^XY'^Z^^ - 288P^v'^w^{80Xiy^ - 3a^)Y^Z^ - 6aP^w{280Xiy^ - a'^)Y'^Z'^ 

+ 15a/3^(40Xi/2 _ a^)XYZ^ - 1882A0v^Paw"^XYZ"^ - 2160v'^P^a^wXYZ^ 

- 6UAiyW{820Xiy^ + a^)YZ‘^ + 5P^{160Xi^^ - a^){20Xv'^ + cP)XZ‘^ 

+ 80u2u,2(128Xu2 + a^){l60Xiy^ - a'^)Z‘^ + Pa{160Xv^ - a^){560Xu^ + a^)wZ^ 

+ 24iy^P^wY^Z - AOP^iy^XY^Z + AOP^iy^awY'^Z + SBSAiy^P^w^Y'^Z + 98304ZY^iy^w^ 

+ 1536iy'^Paw^Y^Z-24iy^P^w(320Xiy^-a^}Y^Z-60iy^aP^XY^Z-8Paiy^w(160Xiy^-a^)Y^Z 
+ 41/2(160X^2 _ a2)(5/32x - 256iy^w^)Y^Z + 4096iy^w^Y^ + 128iy‘^aPY^w + WOXiy^P^Y^, 

where we have written P = u — 4 and a = v — 2. 

Proof. We specialize Theorem 6 to g = 4. The polynomial D{w,x) becomes a square, and is 
independent of w: 

D{w, x) = {2vx + /3)2. 

We we use extensively the notation P = v— 1, a = v — 2. Let us specialize to g = 4 the five 
non-differential identities (88-92). Equation (88) holds verbatim. The identities (89) and (90) 
both specialize to 

2vRq = pPi, 

so that R{w^ x) is a multiple of y/~D(w~v) = 2vx -I- p. However, we can also derive a second 
relation from (89) and (90): if we first eliminate Rq between them, and then set q = 4, we obtain 

{16v^w + aP)Ri - 4vPQi + Si^^Qo + 4vP^ = 0. (107) 

We also have the two relations (91) and (92) relating Ti to the other unknown series. 

We now go back to the differential system (70). We will expand its numerator in powers of 
X = X + P/{2v). We write accordingly 

P{w,x) = X^ + P 2 X‘^ + PiX + P^, Q{w,x) = 2vX'^ pQiX+ Qo, R{w,x) = RiX, 


where the Pj’s and Qj are series in w with rational coefficients in v. Note that we have six 
unknown series instead of seven because R{w,x) is a multiple of X (in other words, Rq = 0). 
The remaining four generic identities, namely (88), (107), (91) and (92), translate as follows: 


Qi = vP2 - 3a/4, 

8iPQq= — {l&v^w + aP) Ri, 

i/Pi = 10z/(2Ti - u>) -k Qo + vP‘^14 Y SaPa/S + ba^/{<o4v), 
Pi = 4u(2r( - 1 ). 


The following initial conditions are translated from (69): 


A(o) = -0. A(o) = M, = 


(108) 

(109) 

( 110 ) 
( 111 ) 


16u3 ’ 


Pa 

2 P 


Qi(0) = -a- p. 


( 112 ) 
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Expanding in X the numerator of (70) gives a system of five differential equations in the six 
series Po,Pi,p 2 ,Qo,Qi,Ri- 

2PoQo - QoPo + 2-Poi?i = 0, 

2Q'oPi + 2PoQ'i - PoQi - PlQo + 3Ai?i = 0, 

2Q'oP 2 + 2PiQ[ - P[Qi - P^Qo - 2vPl, + 4P2i?i = 0, 

2Q'i A - Q 1 P 2 - “i-vP'i + 2Qo + 5Pi = 0, 

Q\ - lyP' = 0. 


We will need to inject an additional equation (namely (109)) to complete this system. 

The fifth (and last) equation of the system is solved by (108). 

The fourth equation is a consequence of (108), (110) and (111). 

In the first equation of the system, let us replace Pi by its expression in terms of Qq derived 
from (109). Using the initial conditions (112), the resulting equation (in Pq and Qq) can now be 
integrated into: 

(I 6 + a/3) Pq = —/SQq. (113) 

Let us finally eliminate Pi between the first and second equations of the system. The resulting 
equation can be integrated into 

PiQo . 


2Qi = 

for some constant c. Using the initial conditions (112), this constant is found to be zero, so that 

2 Q 1 P 0 = PiQo. (114) 


This is as far as we have been able to go in the integration of the system. At this stage, we have 
seven unknown series (the Pj’s, Q^ ’s, Pi and Ti) related by one differential equation (the third 
one in the system), and six “non-differential” equations (namely (108-111), (113) and (114)) in 
which the only derivative is T[. Using in the following order (108), (111), (109), (110), and (113), 
we can now express Qi, Pi, Qoi Pi^ Po terms of P 2 , T[ and Pi. Plugging these expressions 
in the remaining non-differential equation (114) and in the remaining differential equation gives 
a pair of differential equations between P 2 and Pi. Given the form of (110) and (111), it makes 
sense to write them in terms of S' = 2Pi — w. They read: 

4 (AI3vp2 + + 5a/3) S' = IGz/^Pf 2Avap2 + + 5a^ 


2 + ap) (Avp 2 -f 3a - 2/3S') S" - 32z/2/3S'^ 

-f 81 / ((iGi/^iu -H a/3) P 2 -k Av‘^p 2 - S' - {l2v'^ap2 + 3201/^5 -k P^ = 0. 

We finally eliminate P 2 as follows: we differentiate the first equation above to obtain a total of 
three polynomial equations in S,S',S", P 2 and P 2 , from which we eliminate P 2 and P 2 using 
resultants. This yields the differential equation of the theorem. 

To finish, one checks that the recurrence relation that underlies the differential equation for 
Pi determines all coefficients of this series once we prescribe the first three coefficients of Pi. g 

8.2. General planar maps 

Theorem 14. Let Mi = Mi{h',w,t) be the generating function of four-coloured planar maps, 
counted by vertices (w), edges (t) and monochromatic edges (p) (as before, the root-vertex is 
coloured in a prescribed colour). This series is characterized by the initial conditions Mi = 
w -k w(v -k iw -k vwp -k 0{P) and a differential equation of order 3 and degree 11. 

The differential equation can be seen in the Maple session accompanying this paper, on the 
authors’ web pages. 
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Proof. We specialize Theorem 1 to g = 4. The polynomial D{t,x) is again a perfect square: 

D(t, x) = ((u + l)x — 2)^ . 

Equation (94) still holds of course. The identities (95) and (96) both specialize to 

{v + 1) i?o + 2Ri + 4(1 - 2w) = 0. 

Given that i ?2 = (u + 1)(1 — 2w) (this is one of the initial conditions in Theorem 1), this means 
that R{t,x) has a factor \/D = (u + l)x — 2. As in the case of triangulations, we can derive a 
second identity from (95) and (96) by eliminating Ri between them: 


~ (i® T 2) (/3 + dm) ifo + (,5 + 2) (/3 + dm) Qq + 2 (/3 + Aw) Qi — 

A{2w-l)ij3 + 2){p + Aw)t-8w, (115) 


with j3 = V — 1. The two characterizations (97) and (98) of Mi = t^Mi still hold. 

We now go back to the differential system (47). We will expand its numerator in powers of 
X = x — 2/{y +1) instead of x. We write accordingly 

P(t, x) = X'^ + P^X^ + p 2 X^ + PiX + pQ, Q{t, x) = X^ + QiX + Qq, 


and 

R(t, x) = (iy + 1)(1 - 2w)X^ + RiX, 

where the Pj, Qj and Rj are series in t with rational coefficients in u and w. We have taken into 
account the fact that i? is a multiple of X. The remaining four identities, namely (94), (115), 
(97) and (98), translate as follows: 


Qi =P3/2 + 2(2u;-l)(/3 + 2)t, 


2((/3 + 2)^(/3 + 4w)t — /3)i?i + (/? + 2)^(/? + 4w)Qo — 0, 
p 2 = 12(13 + 2)"^Ml + 2Qo + Pi/A - 2(2w -!)(/? + 2 )tp3 - 12tw/3 - I2t, 
Ri = 2(13 + 2)‘^M[ - {2w - l)(/3 + 2 )p3/2 - 2/3w - 2. 


Eq. (118) actually results from the elimination of Qi between (97) and (116). 
The following initial conditions are translated from (46): 


PoiO) = 


4/32 

(2 + ^)4’ 


Pi(0) = 


W - 2 ) 

(2 + /3)3 ’ 


Qo(0) 


2/3 

(2 + / 3 ) 2 ’ 


Qi(o) 


2-/3 

2 + 13' 


(116) 

(117) 

(118) 
(119) 


( 120 ) 


Expanding in X the numerator of (47) gives a system of (only) six equations in the seven 
unknown series Pq, Pi, P 2 , P 3 , Qo, Qi, Pi- 


2P0Q0 - QoPq + 2P0-R1 = 0 , 
2Q'oPi + 2PoQ'i - PoQi - P[Qo + 3 PiPi = 0, 
2 Q'op2 + 2 PiQi - P[Qi - P2Q0 -Po + 4 P 2 i?i - ( 2 w - l ){(3 + 2 )Pi = 0 , 
2 Q'iip 3 + 2p2Q'i - P2Q1 - P3Q0 -P[+ 5P3R1 - 2 { 2 w -l){f 3 + 2)p2 = 0 , 
2Q'iP3 - Q1P3 -P2 + 2Qo + 6P1 - 3 ( 2 u; - 1)(^ + 2)p3 = 0 , 

2Q[ - P( - A(2w - l)(/3 + 2) = 0. 


We will need another equation (namely (117)) to complete this system. 

The sixth (and last) equation of the system is solved by (116). 

The fifth equation is a consequence of (116), (118) and (119). 

In the first equation of the system, we now replace Pi by its expression in terms of Qq derived 
from (117). The resulting equation (in Pq and Qo) can now be integrated into: 

((/3 + 2)2(/3 + Aw)t — [3 )Pq + /3Qo = 0. 


We have used the initial conditions (120) to determine the integration constant. 


(121) 
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Let US finally eliminate i?i between the first and second equations of the system. The resulting 
equation can be integrated into 

on — I /l^ 

2Qi — —^-h c\J Pq, 

Pq 

for some constant c. Using the initial conditions (120), this constant is found to be zero, and 

2QiPo = PiQo- ( 122 ) 


This is as far as we have been able to go in the integration of the system. Using (in this 
order) (116), (119), (117), (118), (121), and (122), we can now express Qi, Ri, Qo, P 2 , Pq and 
finally Pi in terms of P 3 , M[ and Mi. Plugging these expressions in the remaining equations 
of the system (the third and fourth) gives a pair of differential equations between P 3 and Mi. 
Both have order 1 in P 3 and order 2 in Mi. Denoting A = 2 + /3 and p = A^f(/3 + 4w), they read: 


16A^ {p - /?) (tM'i - 3Mi) M" - pP|M" + 4Xt{2w - 1) {p - 3/3) P 3 M 1 

+ m {tX‘^{2p‘^w + P{2w + if + 8w) - 2/3(/3w + 1)) + SA^p {M[f 

- 4tX{2w - 1) (2p - 3/3) P^Pp + p (#3 - 2A(2u> - 1) j + 12A(2u; - 1) (p - P) P 3 M 1 

- StA^ (4/3 u;2 + 2w{P‘^ + 2/3 + 4) + /3) - twA^ (^Pj + 4At(2w - 1)) P( 

+ 2twX^{2w — 1 )P 3 + 8 twX^{l + w/3) = 0, 

and 

(pPg/A - Ui 2 w - l)(p - 2/3)) /Uf 

+ p (Pg/A + 12 u; - 6 ) M( - 3p/t (Pg/A + dw - 2 ) Mi + 2 twXp 3 = 0 . 

By eliminating P 3 between them, we obtain a differential equation of order 3 for Mi (we first 
eliminate Pg between the two equations to obtain one equation between P 3 , Mi, M[ and M"; 
then we follow the steps described at the end of Section 8.1). 

To finish, one checks that the recurrence relation that underlies the differential equation for 
Ml determines all coefficients of this series once the coefficients of P and P in Mi are prescribed. 


9. Spanning forests of cubic maps: q = 0 

When <7 = 0, the series Ti has an interesting combinatorial interpretation in terms of spanning 
forests, as discussed in Section 2.3. Moreover, the polynomial D{w, x) involved in the differential 
system for triangulations simplifies greatly, and this will allow us to characterize Ti by a second 
order equation. We say a map is near-cubic if all non-root vertices have degree 3. 

Theorem 15. Let G{j3,w) he the generating function of near-cubic planar maps with a root- 
vertex of degree 1, equipped with a spanning forest, counted by the number of faces (w) and the 
number of connected components of the forest, minus one (fd). This series is characterized by 
the initial conditions G = 0{wp and the following differential equation, written in terms of 
W = 2G- w/P: 

0 = - P^{5WP -Pw-1 - w)W'^ + 4{P + l){5WP- Pw + W" 

- 48 P^w(P -h 1)IU'^ -k 8/3(/3 -k 1)(5 WP - Pw-h w)W'^ +4(P^- 1)(5 WP - Pw-I- w)W'. 

Remark. Theorem 15 has been rederived in [10] using a simpler, and more combinatorial, 
argument. 
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Proof. Let us apply ( 8 ) to the set Ai of near-triangulations with outer degree 1. We obtain: 

= (123) 

where /3 = u — 1 and G{j3, w) is defined above. 

Hence let us specialize Theorem 6 to g = 0. The polynomial D{w,x) has degree 1 in x, and 
is independent of w: 

D{w, x) = /3^(1 -I- 4a;). 

Let us now specialize the five generic non-differential equations (88-92). Equation ( 88 ) holds 
verbatim. The identities (89) and (90) give i?o = —2;d and Ri = —8/3, so that: 

R{w,x) = -2I3{1 + Ax) = -2D{w,x)//3. (124) 

Equation (91) does not involve Ti, but it expresses Pi in terms of Qo and Qi: 

- 48u/3w; - 4 Piu 2 -h 4uQo + (Qi - l)(Qi + u - 3) = 0. (125) 


Finally, (92) gives again Pi = —8/3, and we still need to relate Ti (or equivalently T 2 = Ti/v) 
to the series Pj and Qj. This can be done by specializing (93) to g = 0: 

8/3Ti7i/- 4 P 2 -h 1 = 0. (126) 

We now go back to the differential system (70). We will expand its numerator in powers of 
X = a; -I- 1/4 instead of x. We write accordingly 

P{w,x) = X^ + P 2 X^+PiX + Po, Q{w,x) = 2iyX'^ + QiX + Qo, R{w,x) =-8f3X, 

where the Pj’s and Qj are series in w with rational coefficients in u. Note that we only have five 
unknown series instead of seven, because P is completely determined (see (124)). The remaining 
three generic identities, namely (88), (125) and (126), give 

Qi = UP 2 + 1 - u/2, (127) 

4j/Qo = 47Pi -k 27P2 - Qi - 4/3Qi -k A8i^j3w - (7v - 6)(u - 2)/4, (128) 

and 

2vP2 = 4:PT[ - V. (129) 

We will need the values of several series at in = 0, which follow from the initial conditions (69): 

Po(0)=0, Pi(0) = l/16, P2(0) = -1/2, Qo(0) = (u-2)/8. (130) 

Expanding in X the numerator of (70) gives a system of five differential equations between the 
five series Pq, Pi, P 2 , Qo and Qi: 

2 P 0 Q 0 - QoPo = 0, 

2Q'oPi + 2 PoQ'i - -PoQi - HQo - 8/3Pi = 0, 

2 Q'oP2 + 2PiQ'i - P(Qi - P^Qo - 2uP/ - I 6 / 3 P 2 = 0, 

2Q[p2 - QiP^ - 2vP[ + 2 Q'q - 24/3 = 0, 

Q\ - up' = 0. 


The fifth (and last) equation is solved by (127). 

The fourth one is a consequence of (128) and (127). 

The first one is readily integrated into Pq = cQq, for some constant c. The initial condi¬ 
tions (130) give c = 0: in other words, Pq = 0 and P{w,x) is a multiple of D(w,x). 

In the third equation, let us replace Pq by 0, then Qo by its expression derived from (128), 
and Qi by its expression (127). Finally, let us introduce the primitive P of Pj: 


/ 


P2 dm = 



u 


w 

2 


S := 


(131) 
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(the second equality follows from (129), with Ti(0) = 0). The resulting equation (in Pi and S) 
can now be integrated into: 

- 2 (21/5" + iy-2)Pi+ lyS'^ + (i^ - 2) /2 + 48w(3S' - 80/35 = 0. (132) 

We have used the initial conditions (130) to determine the integration constant. 

This is as far as we have been able to integrate the system. We can now express Pi, P 2 , Qi, 
Qo in terms of S and S' using (in this order) (132), (131), (127), and (128). We plug these 
values in the second equation of the system, together with Pq = 0, and this gives us a second 
order equation for S, or, equivalently, for Ti: 

p(ri,r(,T") = o 

with 

P{X, Y, Z) = + 8l3^w{l3 - 8)ZY^ - 80(3^XZY^ - I2l3'^w{p - 2)ZY‘^ + UO^^XZY'^ 

+ QP^wZY - mp^XZY + AQQP^vZX'^ + lO^ZX - 8Ou0^w{v + Z)ZX + AI3{v + ifvw^Z 

- vwZ - l<d2vwl3'^Y^ - lQv(3'^w{l3 - 14 )^^ + imp'^vXY'^ 

- Avj5w{P‘^ -13+ 16)r + 40^/32(/3 - 5)XY + 2v^wl3 - 20^/3(/3 - 3)X. 

We now return to (123) to derive an equation for G{f3, w), or equivalently for W = 2G — w/[3. 

To finish, one checks that the recurrence relation that underlies the differential equation for 
G determines all coefficients of this series once we prescribe G = 0{w^). g 

10. The self-dual Potts model on general planar maps 

To finish, we derive a third order differential equation for the self-dual Potts model defined 
at the end of Section 2.3. 

Theorem 16. Let S = S{/3, t) be the generating function for the self-dual Potts model on planar 
maps. More precisely, let 

5 (^,^)=/3^2M(/3^/3-fl,/3-^^;l) = ^^"(^)+2 ^ ^2c(S)+e(S)-i-v(M)^ 

M SCE(M) 

where the first sum runs over all rooted planar maps M. This series is characterized by the 
initial conditions S = P + 0{P) and the following differential equation of order 3 and degree 4; 

P{S, S', S", S'") = 0, 

where 

P{X, Y, Z, T) = 2p{Xj3Y - 2)(2/3y - 1)(3A/3X - jdt- Af)T - /3V(3A/3A: - jdt- At)Z^ 

- pP{-2Xp^pY^ + /3(16At -P-6)Y + 2AX‘^l3^XY + 6A/3(/3 - 10)X - 2A/3t -f 48t -f 2)Z‘^ 

+ 2{Xj3Y - 2){2Xl3^pY‘^ - P{8ph + 160/3t + l3 + 288t - 10)y -f A8X‘^P^XY - A8Xj3^X 
+ 8ph + 80/3t + P + 6At- A)Z - 8(2/3^ - l)(Ay - l){Xj3Y - 2)^, 
with X = /3 + 2 and p = 8Xt — 1. 

Proof. We specialize Theorem 1 to g = /3^, w = 1//3. The argument of Section 4.6 proves that in 
this case as well, the differential system and its initial conditions define uniquely the series Pj, 
Qj and Rj. Observe that one of the initial conditions reads R 2 = 0, so that R{t, x) has degree 1 
in X. This is the key simplification in the self-dual Potts model. 

The polynomial D{t,x) does not simplify drastically. With the notation of the theorem, it 
reads: 

D{x, t) = A/3^x2 + {I3- 2)pI3^/A, 
with X = x — A/2. We write accordingly 

P{t, x) = P{t, X)=X‘^ + PsX^ + P 2 X‘^ + PiX + Po, 


(133) 
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Q{t,x) = Q{t,X) = + QiX+ Qo, and R{t,x) = R(t, X) = RiX + Rq. (134) 

The initial conditions at t = 0 are even functions of X: 

P(0, x) = {X^ - 1/4)2, x) = X^- 1/4. 

The other two initial conditions read P 4 = 1 and R 2 = 0, and we have already taken them into 
account in (133) and (134). But then it is easy to check that (P(t, —^), Q{t, —X), —R{t, —X)) is 
another solution of the differential system, satisfying the same initial conditions. The uniqueness 
of the solution implies that P and Q are even functions of X, while R must be odd. Hence 

Pi = P 3 = Qi = Rq = 0, 

and we are left with four unknown series only. 

With these values, the generic equations (94) and (96) automatically hold, and the other three 
reduce to: 

U\Qo+piRi-13 + 2) =0, 

< UPXS + 2Qo -P 2 - 4(/l + 4)t = 0 , (135) 

[-2/3A5" + Pi+4 =0, 

where S = (3Mi = pt^Mi is the series defined in the theorem. 

Expanding in X the numerator of (47) gives a system of four equations in the four series Pq, 
P 2 , Qo, Ri, corresponding to the coefficients of X^^ for i from 0 to 3. The other coefficients 
vanish due to the parity properties. 

In these four equations, let us replace P 2 , Qo and Pi by their expressions in terms of S 
derived from (135). The fourth equation then vanishes. The three others are found to be 
linearly dependent (over Q(;5)), and we consider those corresponding to the coefficients of 
and X'^-. 

4(p/3P" + AXPS' - 8)Po = p{2pS' - 1)P/ 

and 

4P/ = ( 2 / 32 / 5 " _ p^(48A^S' - m - 48t - 1))S" + 8XI3^pS'^ - 4/3{P + 6 )pS' + + 128t - 8 . 

It remains to eliminate Pq to obtain the differential equation of the theorem. To conclude, one 
checks that the underlying recurrence relation on the coefficients of S{t) determines all of them 
once the coefficient of is prescribed. ■ 

11. Final comments 

At the moment, our solution does not solve the important question of locating phase tran¬ 
sitions and finding critical exponents of the Potts model on planar maps. Partial results are 
known, for instance for the Ising model [ 8 , 32, 2], for maps equipped with a spanning forest [2], 
and for properly coloured triangulations [35]. The latter work exploits Tutte’s differential equa¬ 
tion (1). Moreover, a parametrized description of the critical value of v (depending on q) for 
planar maps is given in [ 6 ], for 0 < g < 4. The analysis builds on the catalytic variable that 
records the degree of the outer face rather than on the size variable, but this should not influence 
the location of the critical point. Also the results of [23] involve the dependence of the series in 
this catalytic variable. 

It is natural to ask if the Potts generating function might also satisfy a linear differential 
equation. This would make the analysis of critical points much easier, since finding the location 
and nature of the singularities would then be (almost) automatic [24]. However, this possibility 
has been ruled out (for triangulations) by an asymptotic argument in [ 10 ] . 

The second order equations that we have obtained for triangulations in special cases (Sec¬ 
tions 7 to 9) probably deserve a further treatment. Do they fit in Painleve’s classification? 

Finally, our proofs are very heavy, and it would be nice to establish differential algebraicity 
by simpler, and ideally more combinatorial means. This has been done in several special cases 
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(Ising [13, 14], spanning trees [34], spanning forests [10], bipolar orientations [5, 9]). Two par¬ 
ticularly attractive problems are three-coloured maps (known to be algebraic [ 2 ]) and properly 
coloured triangulations (Tutte’s recurrence relation (1)). 
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